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ABSTRACT 


The development of new high-speed, high-capacity public air transportation is 
one of the challenges facing commercial aviation today. One possible solution to this 
challenge is a supersonic oblique flying wing aircraft. This project involves the 
design and construction of a radio-controlled model of a flying wing aircraft as a 
means of examining some of the aircraft’s unique problems, especially its natural 
instability. The different areas of stability, controls, aerodynamics, structures, and 
propulsion are analyzed in the context of the aircraft as a whole. The design relies 
heavily upon computer modeling and analysis to overcome the technical challenges 
encountered. Since the size of the model prohibits the use of automatic control 
systems, the wing is designed to be as inherently stable as possible, and an airfoil is 
designed to provide the necessary characteristics. Finally, the optimum design is 
chosen and constructed using composite materials. 
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EXECUTIVE SUMMARY 


The main challenge of this project was to design an aircraft that will achieve 
stability while flying without a horizontal tail. The project focused on both the 
design, analysis and construction of a remotely piloted, elliptical shaped flying wing. 

The design team was composed of four sub-groups each of which dealt with 
the different aspects of the design, namely aerodynamics, stability and control, 
propulsion and structures. Each member of the team initially researched the 
background information pertaining to specific facets of the project. Since previous 
work on this topic was limited, most of the focus of the project was directed towards 
developing an understanding of the natural instability of the aircraft. 

Once the design team entered the conceptual stage of the project, a series of 
compromises had to be made to satisfy the unique requirements of each sub-group. 

As a result of the numerous calculations and iterations necessary, computers were 
utilized extensively. In order to visualize the design and layout of the wing, engines 
and control surfaces, a solid modeling package was used to evaluate optimum design 
placements. 

When the design was finalized, construction began with the help of all the 
members of the project team. The nature of the carbon composite construction 
process demanded long hours of manual labor. The assembly of the engine systems 
also required precision hand work. 

The final product of this project is the Elang, a one-of-a-kind remotely piloted 
aircraft of composite construction powered by two ducted fan engines. 
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aspect ratio 
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chord length 
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section drag coefficient 
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total drag coefficient 
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location of the aerodynamic center along the y-axis 
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location of the center of gravity of added mass along the y- 

gr.2 


axis 


Yweg 

location of the center of gravity of the wing along the y- 

— 


axis 


a 
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CHAPTER ONE: INTRODUCTION 


The development of new high-speed, high-capacity public air transportation is one 
of the challenges facing commercial aviation as our society enters the twenty-first 
century. One possible solution to this challenge is a supersonic oblique flying wing 
aircraft. For instance, such an aircraft could carry as many as 500 passengers travelling 
at 1100 miles per hour, enabling transportation from New York to Los Angeles to be 
completed in under three hours (Jones, "Technical Note", p. 106). This concept was 
initially proposed by Robert T. Jones in 1958 (p. 104) but was not believed to be a 
feasible design due to structural and stability considerations. Today, however, the 
technology exists to overcome these challenges. 

Since the full scope of such an endeavor involves numerous difficulties, simplicity 
must be used to surmount these problems by examining the different aspects individually. 
With this approach in mind, a simplified design was developed and analyzed which 
studied the stability, structural and aerodynamic issues of an unswept elliptically shaped 
flying wing. Even the nature of this type of aircraft provides design obstacles which 
must be faced in order to eventually apply the solutions to a full scale aircraft. 

The result of these efforts is the Elang. A Figure 1.1 (also Appendix G) shows 
the final design. Due to the unique design requirements, the development of an original 
airfoil was necessary, resulting in the HPRS 33 airfoil section. The wing has a root 
chord length of 45.5 cm and a span of approximately 2.5 m, with an aspect ratio of 7. 
The planform area is elliptically shaped with a straight quarter-chord line according to, 
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c(z) = = °- 455 / 1 -( T^ )2 

The induced drag is minimized by using an elliptical shape giving a total drag coefficient 
of 0.02208. The flying wing is at a geometric attack angle of 4 degrees and as a results 
generates a lift coefficient of 0.4355 at cruise speed of 15 m/s. The lift to drag ratio is 
therefore 19.69. 

The main control surfaces include two ailerons and one elevator which are used 
to stabilize the aircraft in flight. The ailerons are located on each end outboard the 
engines, occupying 2.9 percent of the wing planform area. The elevator area is 13 
percent of the wing planform area. 

Two vertical tails are utilized, each located on the rear portion of the engine 
ducts. The airfoil section is the symmetric SD 8020 airfoil. The vertical stabilizers each 
have an effective area of 0.10716 m 2 to provide yaw control. Rather large rudder 
surfaces were purposely employed to ensure adequate control. Only vertical tails are 
used, since by nature the flying wing design does not contain a horizontal tail. 

A propulsion system was chosen which provides ample thrust and relatively clean 
aerodynamics. The system consists of two ducted fan engines, each placed 57 cm from 
the root chord. Operating at 22000 rpm, each engine provides 15.60 N of thrust. One 
centrally located fuel tank supplies both engines. 

The wing structure is designed to be composed of a foam core wrapped with a 
skin of carbon fiber and epoxy. This provides a relatively lightweight structure that can 
still withstand sizable loads. 
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The main focus of this project involves the design and construction of a radio 
controlled model of an elliptical flying wing as a means of examining some of the 
aircraft’s unique problems, especially its natural instability. The design relies heavily 
upon computer modeling and analysis to overcome the technical challenges encountered. 
Since the size of the model prohibits the use of automatic control systems, the wing is 
designed to be as inherently stable as possible, and an airfoil is designed to provide the 
necessary characteristics. Finally, the optimum design is chosen and constructed using 
composite materials. To achieve this final product the group was divided into four sub- 
groups, dealing with aerodynamics, controls and stability, structures and propulsion. In 
the following chapters, the analysis of each sub-group and explanations for its results are 
presented. 



Figure 1.1: The Elang 
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CHAPTER TWO: STABILITY AND CONTROL 

SURFACES 

2.1 Stability 


Stability is defined as the ability of an airplane to fly straight with its wings level, 
and remain in equilibrium without pilot assistance. It can be divided into two categories, 
static and dynamic, and can be analyzed along three axes, pitch, roll, and yaw. 

Static stability deals with the natural ability of an aircraft to create a restoring 
moment after a disturbance, and return towards its initial equilibrium position. For 
instance, when in straight and level flight a statically stable airplane encounters a gust 
of wind, it generates a moment to counteract the perturbance. This type of stability 
does not give any information on how long it takes to return to its original position, and 
does not even insure that the plane will return exactly to its equilibrium state. A 
dynamic stability analysis must be performed in order to obtain this information. 

Dynamic stability concerns the motion of the aircraft after the disturbance. In 
order for the aircraft to be dynamically stable, it must return to the equilibrium position 
within a certain amount of time, and thus damping is required. If an aircraft is 
dynamically stable, it will also be statically stable, but the reverse is not necessarily true. 
There are three axes about which stability must be achieved. The rotations are 
about the conventional x, y, and z axes, and the motions are designated as roll, pitch, 
and yaw, respectively (Figure 2.1). These motions can be separated into two distinct 
groups, lateral and longitudinal. Lateral motions consist of rolling, yawing, and 
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sideslipping. Longitudinal motions denote the 
movement along the aircraft’s flight path, its 
vertical movement, and pitching. 

The major difficulty in stabilizing a flying 
wing is along the pitch axis. In conventional 
designs, a horizontal tail offsets the moments 
generated by the wing and fuselage. However, a flying wing does not have a horizontal 
tail. Consequently, an alternative, more advanced method for stabilizing the wing had 
to be utilized. 



Figure 2.1: Axes of Rotation 
(Shevell, Fundamentals of Flight) 


2.2 Stability Analysis 

Throughout the analysis of stability and controls, the major tool that was utilized 
was the MATHCAD software package. By using this package, changes could be made 
to the iterative analysis easily and accurately. A listing of the MATHCAD results is 
located in Appendix A. 

2.1.1 Pitch Stability 

Conventional airfoils typically have a negative moment coefficient about the 
aerodynamic center. After analyzing a few such airfoils, it was found that to achieve 
static stability, the engines would have to be placed about a meter below the wing with 
some airfoils, and at the back of the wing for others. Since the back of the airfoil is the 
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thinnest part, and one meter is considerably too long for clearance, neither of these 
designs was acceptable. Thus, it was determined that an airfoil with a positive 
coefficient of moment about the aerodynamic center was needed. The positive moment 
causes the wing to have a natural tendency to pitch up. Therefore, the aerodynamics 
sub-group designed a new airfoil, through an extensive iterative process, which has a 
high lift to drag ratio along with a positive moment about the aerodynamic center. 

Once the airfoil was selected, engine placement was used to obtain stability. The 
center of gravity for the design was calculated from the known values of the components. 
The weight of the wing and the location of its center of gravity along the x and y axes 
were generated using ARIES, a solid modeling computer program. These values were 
found to be 16 N, 0.209 m, and 0.00561 m respectively . The engines each weighed 
5.719 N. A fixed mass of 0.5 kg was added to allow for the receiver, batteries, and 
other components. This mass was placed 5 cm from the leading edge along the chord 
line. The center of gravity of the design was calculated by Equations 2.1 and 2.2. 


x (2 ) +* *w +jc *w 

_ *ct e W„ w mass mass 


eg 


2*w +>v 

e w mass 


[ 2 . 1 ] 


V = « 

J cg 




2*w> +w +w 

e w mass 


[ 2 . 2 ] 


The moments about the center of gravity of the wing were analyzed, and the 
results led to Equation 2.3. 
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An iterative process was conducted with Equations 2.1 and 2.3, to determine the 
placement of the engines along the x-axis to ensure no inherent pitching moment. 
Equation 2.2 was used to find the location of the center of gravity of the wing along the 
y axis. Initially the y location of the engines was estimated to be six millimeters to begin 
the iteration of the design. The calculations were carried out with lift, drag, and moment 
data for the cruising angle of attack of 4 degrees. 

A Moment Coefficient vs. Angle 
of Attack graph was generated to 
determine whether or not the wing was 
statically stable (Figure 2.2). By 
definition, an aircraft is statically stable if 
it satisfies Equation 2.4. 



Figure 2.2: C M vs. a 


dC 


M 


da 


<0 


[2.4] 


As seen in Figure 2.2, the negative slope ensures static stability along the pitch axis. 

A major concern with this flying wing is the very low moment of inertia about 
the pitch axis, which is a result of the wing having no sweep angle. This low moment 
of inertia indicates that there is little time for a pilot to react to disturbances during 
flight. This problem, coupled with inability to guarantee exact placements and weights 
of the internal components of the wing, makes controllability extremely difficult for the 
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pilot. Therefore, a gyroscope system was included as a safety precaution. The 
gyroscope monitors deviations about the pitch axis, and makes constant adjustments by 
using the elevator. The gyroscope can be utilized at the pilot’s discretion, depending on 
flight conditions. 

In reference to dynamic stability, all available documentation depended on a 
horizontal tail to provide the damping. However, the flying wing has no horizontal tail. 
Thus, any damping that would occur to provide dynamic stability along the pitch axis is 
quite small. The lack of damping again showed the need for the gyroscope system. 

2.2.2 Roll and Yaw Stability 

The stabilities in the roll and yaw directions are linked closely together. Rotation 
along one of these axes induces an accompanying moment along the other axis. Due to 
the symmetry along the x and z axes, the problem of stability in these directions is not 
anticipated to be present. Additionally, if the wing does roll or yaw, the moment of 
inertia about these axes is substantial enough to create a rate of rotation that will provide 
ample time for the pilot to counter the deviations. 

2.3 Control Surfaces 

The control surfaces consist of an elevator, a pair of ailerons, and two rudders, 
which control motion along the pitch, roll, and yaw axes, respectively. The analysis was 
performed by considering the lift, drag, and moment characteristics of each control 
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surface. Each control system was analyzed independently of the other two. When sizing 
the control surfaces, they were purposely overdesigned, because the controllability of 
previous design projects has proven to be inadequate. 


n 


It 



2.3.1 Elevator Analysis 

The data for the elevator was obtained by modifying the inviscid code to provide 
new lift, drag, and moment coefficients when the control surface was deflected. The 
analysis was carried out through a range of angles from -30 to 30 degrees. 

Once the data for the deflected area of the wing was procured, it was combined 
with the previous data for the remaining portion of the wing. From the combination of 
this data, a modified moment equation was developed. For these calculations, the length 
of the elevator was 1.4 meters at 18 percent of the cord length, resulting in an elevator 
area equal to 13 percent of the wing planform area. The moment is given as: 


-yJ * D A S ~ D ^ -yJ 

Vi -*« Vv "W * D ’<&* " W - W 

where subscript def denote characteristics 
from the deflected portion. 

Once the moments were 
determined, Newton’s Second Law was 
applied to calculate the pitching 
acceleration of the wing. Figure 2.3 



Figure 2.3: Angular Acceleration vs. 
Elevator Deflection Angle 
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shows the results of this analysis. The values for this acceleration ranged from negative 
to positive 200 rad/sec 2 . These large accelerations develop because of the small moment 
of inertia of the wing about the pitch axis. From this data, it was decided that the 
gyroscope system should be utilized to aid the pilot in controlling the wing. 


2.3.2 Aileron Analysis 

The data from the elevator analysis was utilized in the analysis of the ailerons. 
This similarity occurs because the elevator and the ailerons deflect along the same axis, 
and both occupy 18 percent of the cord length. When examining the ailerons, it was 
decided that they would deflect at the same rate in opposite directions. For example, if 
one aileron was deflected up 10 degrees, the other one would deflect down 10 degrees. 
The rolling moment created by the drag from the rudders was neglected. The analysis 
produced the following moment equation. 




[ 2 . 6 ] 


Through iteration it was determined that the each aileron would have a length of 
21.8 cm, which resulted in each aileron having 2.9 percent of the wing planform area. 
As with the elevator analysis, Newton’s Second Law was applied to generate the angular 


acceleration about the roll axis of the wing. The acceleration ranged from negative to 
positive ten radians per second squared. 
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2.3.3 Vertical Stabilizer Design 

The most important consideration in the design of the vertical stabilizers was the 
sizing. Due to the inherent control problems involved with a flying wing, it was 
imperative that the control surfaces be oversized to compensate for unpredictable stability 
failures. For most airplanes with a fuselage, tail and vertical rudder, the area required 
for the vertical tail is 3 to 5 percent of the entire planform area of the plane. The design 
of the flying wing called for an area of more than the typical percentage. Thus it was 
determined that the vertical tail area be 12 percent of the wing planform area. With a 
planform area of 0.893 m 2 , the total vertical stabilizer area is 0.10716 m 2 . 

2.3.3 Rudder Analysis 

The inviscid analysis that was used for the elevator and ailerons was also used to 
examine the rudders. Since there is such a small moment arm for the rudders, the size 
of each one is considerably larger than that of conventional designs. Two rudders were 
selected, for symmetry, and for structural considerations. Each rudder will have an area 
equal to 31.7 percent of the vertical stabilizer planform. This sizing is accomplished by 
having the height 21.3 cm, the root cord 17.2 cm, and the tip cord 8 cm. 

In the analysis of the rudder, it was determined that each rudder would be 
deflected the same amount. The moment equations were developed separately for each 
rudder, and then summed. These expressions are given in Equations 2.7 and 2.8. 
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M rr = -D r *z rr -L r (x rr +x tK ) [2.8] 

Once the moment equations were 
developed, Newton’s Second Law was 
applied to calculate the angular 
accelerations. A plot of these 
accelerations is given in Figure 2.4. In 
this graph, the maximum acceleration 
occurs close to zero. Since an inviscid 
analysis was used to obtain the data, drag 
effects are assumed small. In reality however, the drag effects will contribute greatly 
to the moment generated at large deflections. 
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Figure 2.4: Rudder Angular Acceleration 
vs. Angle of Deflection 
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CHAPTER THREE: AERODYNAMIC 
DESIGN AND ANALYSIS 


3.1 Introduction 

This chapter discusses the aerodynamics of the aircraft, including the selection of the 
airfoil, airfoil characteristics, aerodynamic behavior of the aircraft, and the aerodynamic 
shape of the aircraft as a whole. The relationship between stability considerations and 
airfoil design is discussed thoroughly, since stability of the aircraft was of paramount 
importance. The analysis of the lift and drag characteristics for the aircraft turns out to 
be quite simple, since the entire aircraft is an elliptically-shaped wing. 

3.2 Final Design 

The final design is an elliptically shaped wing with a straight quarter chord line, the 
planform view of which is shown in Figure 3.1. The airfoil section used for the wing 
is the HPRS 33, a custom-designed airfoil. This airfoil, presented in Figure 3.2, is a 15- 
percent thick airfoil section with a reflexed camber. The wing has a span of 2.5 meters, 
an aspect ratio of 7, giving a root chord of 0.455 m and planform area S ref of 0.893 m 2 . 



Figure 3.1: Planform of Wing 
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Figure 3.3: SD 8020 Airfoil 


The two identical vertical tails each have a height of 0.213 meter, a root chord of 
0.172 meter, and a tip chord of 0.08 meter. The trailing edge is unswept. The airfoil 
section used for the vertical tails is the SD8020, a 10-percent thick symmetric airfoil. 
This airfoil is shown in Figure 3.3. 

3.3 Aerodynamic Requirements 

The selection of the airfoil and the design of the wing shape were dictated by the 
consideration that the aircraft was to be a flying wing, which had to be controllable in 
flight by a single RC pilot. Three requirements in particular were: the need for 

relatively high lift and low drag; the requisite of static and dynamic stability; and 
structural rigidity as well as the location of engines, servos, and avionics. 

One of the purposes of a flying wing configuration for a transport aircraft is to 
transport people and goods as efficiently as possible. Therefore the aircraft needs to 
have a high lift to drag ratio. A model of such an aircraft is no exception to this 
requirement. Since a finite wing of elliptical planform experiences a minimum induced 
drag, this shape was chosen for the flying wing. The wing quarter chord line is set to 
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be a straight line, visible in Figure 3.1. Having this configuration avoids further 
complications during construction. 

The requirement of static and dynamic stability without a horizontal tail necessitates 
the use of an airfoil with special stability characteristics. While most airfoil sections 
have a negative moment coefficient, the flying wing application requires an airfoil with 
a positive moment coefficient, as discussed in Chapter 2. One way to accomplish this 
is to use an airfoil with a reflexed camber. 

In addition, the wing must withstand the aerodynamic loads, particularly the twisting 
and bending loads. Since the aircraft has no fuselage, all components which include 
receiver and batteries, servos, and fuel tank, must be contained internally in the wing. 
These reasons led to the development of a thick airfoil, namely the HPRS 33 airfoil 
section. 


3.4 Airfoil Design and Selection 


3.4.1 Wing Airfoil Section 

Since the flying wing has such special design requirements, it proved impossible to 
select an existing airfoil with the required lift, thickness, and moment characteristics. 
Therefore, it was necessary to develop an airfoil with the desired properties. To do this, 
the geometries of a thick high-lift airfoil and an airfoil with a positive moment coefficient 
were merged. A two-dimensional inviscid vortex panel code was utilized to calculate the 
lift, drag, and moment coefficients for the resulting airfoil. Points on the airfoil were 
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then changed and entered into the code 
in an iterative process until a 
satisfactory design was produced. The 
end result was the HPRS 33 airfoil, 
shown in Figure 3.2. Figure 3.4 
presents the lift, drag, and moment 
coefficients versus angle of attack. 

Note that the drag coefficient is very 
low, since an inviscid approach was 
used. 

Appendix B lists the airfoil coordinates, as well as the lift, drag, and quarter-chord 
moment coefficients at various angles of attack for the HPRS 33 airfoil section. 

3.4.2 Vertical Stabilizer Airfoil Section 

The first consideration in the design of the vertical stabilizers was the shape of the 
airfoil to be used. After researching vertical tails on conventional airplanes, it was 
determined that a symmetrical airfoil was the best choice. By using a symmetrical 
airfoil, unnecessary moment forces caused by uneven lift distribution could be avoided. 
Research of thin symmetric airfoils provided the SD 8020 shape (see Figure 3.3). 
Appendix C lists the coordinates of the SD 8020 airfoil. 



Figure 3.4: Aerodynamic Coefficients for 
HPRS 33 
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3.4.3 Vertical Stabilizer Sizing 


Sizing is the most important 
consideration in the design of the vertical 
stabilizer or tail as discussed in Chapter 
Two. With a given total vertical tail area 
of 0.10716 m 2 ,the final configuration of 
each vertical stabilizer consists of a root 
chord of 17.2 cm, a height of 21.3 cm, a 
tip chord of 8.0 cm and an aspect ratio of 
0.85 (see Figure 3.5). These control 
surfaces provide adequate area for 

3.5 Aerodynamic Calculations 

The following is a very brief overview of the calculations performed to obtain 
theoretical lift and drag data for the aircraft. 



Figure 3.5: Vertical Tail 


stabilization about the yaw axis. 


3.5.1 Lift 

To obtain theoretical lift values for the aircraft, the slope of the lift curve for the 
airfoil section first had to be obtained. This was obtained simply from the airfoil section 
data given in Appendix B; for the HPRS 33, dc,/da = 6.97 per radian. 


The lift curve slope for the finite wing is then given as 


dC L _ 2nAR [ 3 . 1 ] 

da 2+/(Ai? 2 /tif ) (l+tan 2 A-Af z ) 

The airfoil efficiency. 7\. is approximated as 1.0. The sweep angle, A. is zero, and 
the flight mach number is 0.044. Equation 3.1 thus results in dc L /da = 5 per radian. 
Integrating this expression gives 

C L = 5a + constant [3-2] 

The constant is obtained by using boundary conditions, which are approximated from 
the section lift curve slope: at a = 1°, c, = 0.17. Solving, the constant is found to be 
approximately 0.084. The lift curve slope for the entire wing is 

C L = 5a + 0.084 [3-TJ 

where the angle of attack is in radians. The lift coefficient of the wing can thus be found 
for various angles of attack. 

Knowing the lift coefficient, the lift force can be found at any speed for different 
angles of attack using the equation 


L = 


1 

2 


P v » C^zef 


[3.4] 


The lift of the aircraft at various angles of attack and various flight speeds is 
tabulated in Appendix D, using Equation 3.4. At cruise conditions (4 AoA, 15 m/s), 
the lift coefficient is 0.4355 and the total lift force is 53.55 N. The variation of available 
lift with velocity is shown in Figure 3.6. 
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Because the inviscid vortex panel code was not capable of providing accurate drag 
results, an empirical method suggested by Raymer was used. 

This method depends upon the calculation of the parasitic drag, C Do , based on 
equivalent skin friction coefficients C fc , component form factors FF C , interference effects 
Qc> wetted area S wete , miscellaneous drags C Dmisc , and estimated leakage and protuberance 
drag C D LiP . The parasitic drag is calculated according to: 
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[3.5] 


"'Do 


= E (WgAeJ 


? ref 


+ Co.mlBC + Qj.LP 


(Raymer, pp. 280 -289). Appendix D lists the values of these parameters for each 
component of the aircraft, as well as the total parasitic drag for various flight conditions. 
At cruise conditions (4° AoA, 15 m/s), the approximate value of parasitic drag is C Do = 
0.01346. 


The total induced drag, or drag due to lift, is given as 


C Di = KC , 2 where K = — - — [3.6] 

Dl 1 %ARe 


For an elliptical planform, the Oswald span efficiency factor e is unity. Therefore, 
K = 0.04547, leading to a cruise C Di of 0.008624. The total drag coefficient is simply 
the sum of the parasitic and induced drag coefficients, so that, at cruise, 

C D = Cjte + C Di = 0.01346 + 0.008624 = 0.02208 [ 3 . 7 ] 

The cruise drag force is therefore calculated at 2.72 N and the lift to drag ratio at 


cruise is 53.55/2.72 = 19.69. 



« CHAPTER FOUR: STRUCTURAL DESIGN AND 

ANALYSIS 

E : : -5 

w 


B 


4.1 Introduction 

The structural design of the wing involved two primary objectives. The first was 
to reduce aeroelastic effects since this would alter the distribution of the forces along the 
wing and consequently alter the stability. Therefore the wing was designed to be as rigid 
as possible to minimize any bending. The second was to maintain manufacturability, 
since most of the construction was to be done by students using equipment available on 
campus. 

One can describe the methods of the structures aspect of this project as a simple 
procedure. First, the constants in this procedure were determined. These constants were 
found to be as follows: the material of the skin, a carbon composite weave strengthened 
with an epoxy; the core of the sandwich construction, a polystyrene foam provided and 
cut by an outside vendor; and, finally, the basic interior layout of the aircraft. 

Next, the variables in this procedure which might have caused problems in the 
design were determined. These variables were found to be the values of the stresses, 
forces, etc. that are found through the wing loading analysis. The values which the 
aerodynamics and stability groups produced were substituted into the appropriate 
equations and compared to the safe values of the carbon composite weave. If the stresses 
or shear forces were unacceptable, the other groups were notified in order that the 
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procedure be started over with a new design, until a legitimate final product was 


obtained. 


4.2 Final Design 


Figure 4. 1 shows the final interior design of the wing. The absence of a fuselage 
posed some design limitations, namely placement and access to the electronic and engine 
components. Thus all the components had to be placed within the wing. In order to 
simplify construction most of the internal components were placed together in the mid- 
section of the wing, with access to them through a "trap door" placed on the bottom. 

The engines posed a greater problem since they are the heaviest component on 
the wing and are placed 57 cm from the root creating regions of increased stress 
concentration. These connections were simplified as much as possible to ease 
construction, but still maintain structural integrity. The final engine connection attached 
the duct of the engine directly to the carbon shell by four screws. 

4.3 Method of Construction 

The wing is a sandwich construction with a carbon fiber composite "face sheet" 
and an extruded polystyrene "core". This type of construction makes use of the 
advantages offered by both of the materials used in distributing the load. The "face 
sheet or the carbon composite material will carry the majority of the tension and 
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FIGURE 4.1: INTERIOR DESIGN OF FLYING WING 





compression loading, while the "core" or extruded polystyrene will carry most of the 
shear as well as the compression perpendicular to the skin. 

In order to construct this type of "material", a process called vacuum bagging was 
used since it would give the best results. The principle behind vacuum bagging involves 
placing various layers of carbon/epoxy over the components to be covered. Then an 
evenly distributed force, of magnitude 4-8 Newtons, is applied across the entire surface 
of the part. The result is a reduction in the presence of air bubbles as well as removal 
of excess resin. The presence of excess resin, besides increasing the weight, causes the 
formation of stress concentration points by inconsistency in strength. 

The extruded polystyrene was chosen because it is a rigid foam and therefore 
helped to attain the desired rigidity. In addition the polystyrene is a foam that is easy 
to handle. The foam was initially a single oversized piece that needed to be cut and 
shaped to the desired shape. The airfoil section needed to be hot wired into the planform 
and the taper had to be appropriately shaped using coarse sandpaper. The engine 
connections were made through hot wiring the circular sections at a 2 degree outward 
canter and a 4 degree downward canter. As mentioned previously, the absence of the 
fuselage forced the placement of the components inside the wing, and thus much of the 
mid-section had to be hollowed out. A spar was placed in this area in order help regain 
some of the rigidity that was lost. 

Once the wing had been properly shaped, the carbon fiber skin was placed over 
the foam using the above described vacuum bagging procedure. 
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4.4 Material Selection 


Many important factors had to be considered when selecting the material for the 
wing. The material had to be within tolerance of the necessary yield and ultimate 
strength. It had to be stiff, yet have the ability to resist fracturing due to impact during 
a crash landing. Weight was a main constraint as excessive loads will inhibit the plane’s 
ability to fly. Consequently, the material’s weight and density were also a primary 
concern. The temperature limits also had to be considered in conjunction with the 
engine. If the engines were to become extremely hot, the material had to be capable of 
withstanding the heat so as to resist melting or deforming. Finally, the practicality of 
the material had to be examined. Not only did the material have to be workable with the 
equipment on hand, but cost and availability were also principal factors. 

After much consideration and research, composites were decided to best meet the 
design requirements. Composites met all of the constraints and offered the largest 
reduction in weight. 

A composite is a combination of two chemically distinct phases whose properties 
and structural performance is superior to those of the constituents acting independently. 
A composite is composed of a fiber and a matrix. The fiber supplies the strength and 
stiffness of the composite and carries most of the load while the matrix supports and 
transfers the stresses to the fibers. The matrix also protects the fibers from physical 
damage and the environment. Finally, it provides ductility and toughness to prevent the 
propagation of cracks. Composites cannot take concentrated loads; therefore, the loads 
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should be evenly distributed throughout the composite. Composites are only strong in 
the direction of the fibers; thus, orientation of the fibers is very important. Bidirectional 
fibers at 90 degrees will withstand vertical and horizontal forces while fibers at 45 
degrees are better for torsion. It was determined that the drag forces would be stonger 
than the forces due to torsion; therefore the 90 degree fiber orientation was chosen. 
Carbon was the fiber selected (see Table 4.1). Carbon was selected because of its high 
strength to weight ratio, inherent rigidity, low density, and low cost. Carbon is also 
relatively easy to mold. 


Table 4.1: 

Properties of Reinforcing Fibers 



TENSILE 

ELASTIC 




STRENGTH 

MODULUS 

DENSITY 

RELATIVE 

TYPE 

(MPa) 

(GPa) 

(kg/m 3 ) 

COST 

Boron 

3500 

380 

2600 

Highest 

Carbon 
High strength 

3000 

275 

1 

Low 

High modulus 

2000 

415 


Low 

Glass 

Ml 


MM 


E type 


73 

l 

Lowest 

S type 

— 

85 

B9 

Lowest 

Kevlar 

HM 


MM 

MM 

29 


62 



49 

MM 

117 

WKW 

MM 


Hexel epoxy was chosen as the matrix. Hexel epoxy was chosen primarily for 
its long pot life. This factor made it easier to work with during construction. Epoxy is 
a thermoset plastic meaning that once it is cured it is "set" and the reaction cannot be 
reversed. This made construction a little difficult because no mistakes could be made. 
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Composites are very strong in tension, but they are relatively weak in 
compression. For this reason the composite was supported by foam. The foam provided 
a "cushioning medium" through sandwich construction. Foam did not add much weight 
and was considerably better in compression. Density was the main consideration when 
selecting a foam. The higher the density, the better the compression but the heavier it 
would be. Polystyrene was chosen because it offered an acceptable density and was easy 
to form. 


4.5 Wing Loading Analysis 

One of the most important aspects of the structure of the plane was to 
guarantee minimum bending. If the plane were to bend or twist during flight, the 
distribution of forces may change greatly which would adversely affect the stability of 
the plane; therefore it was extremely important to maintain the shape of the platform 
during flight. In order to verify that the material selected to build the wing would not 
bend under the loading encountered during flight, an analysis of the stresses at the 
wing root was performed. The wing was modeled as a cantilever beam since it is 
symmetric along the root chord. In order to simplify the analysis, the cantilever 
beam was assumed to have a constant area cross-section unlike the actual model 
which has an elliptical taper. 

Since the wing is of sandwich construction, which consists of two materials, it 
had to be "transformed" into a single equivalent material. This transformation was 


4-6 



done by creating a ratio of Young’s Modulus for the two materials involved, in this 
case the extruded polystyrene foam and the carbon fiber composite. This ratio then 
determined an equivalent area for the material. The extruded polystyrene foam was 
converted mathematically into carbon fiber composite. Below is the ratio of the 
Moduli of Elasticity of the two materials. 


n 


J foam 


E, 


car bon fiber 


[4.1] 


Since this ratio is much less than one it was concluded that the foam can be neglected 
in the analysis. 

The loading was then applied to the beam according to those loads estimated 
for level cruising. For example, the lift was elliptically distributed along the span. 
Additionally the drag was uniformly distributed. The weight of the engines and thrust 
were placed in their appropriate position, according to the numbers obtained from the 
stability subgroup. The normal stresses were determined using the following 
equation: 

_ P M x*y A M y* X [4.2] 

AT T 

xx yy 

The results are presented in Appendix E. The maximum stresses encountered at the 
root were calculated to be 23.99 kPa, which is well within the limits the carbon fiber 
composite strength. 
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In addition to normal stress, the 
shear stress was also determined at the 
root. Once again asymmetric bending 
was used. The airfoil section was modeled as a closed loop cross-section with a 
positive shear direction show in Figure 4.2. The results of this analysis are also 
shown in Appendix E. The shear was estimated to be 108.3 kPa, which is once again 
well within the limits of the composite material. 

It was realized that the theoretical material strengths were greater than those of 
the actual values. However, it was believed that the stresses were well within 
tolerance. The objective was to build a rigid, yet manufacturable wing. Using the 
suggested method of construction with the chosen materials, it is felt that this goal has 
been successfully obtained. 



4-8 



CHAPTER FIVE: PROPULSION SYSTEMS 

5.1 Introduction 


Since the flying wing has no fuselage, it was important to find a method of 
propulsion which could be easily integrated with the wing. Many important factors 
needed to be considered when choosing the propulsion system, such as engine type, 
method of propulsion and total system integration. In addition, the propulsion system 
was required to produce adequate and controllable thrust to allow for a smooth flight. 

5.2 Engine Selection 

Two types of engines were considered, electric engines and internal combustion 
engines. Both of these types of engines have advantages and disadvantages with respect 
to their usage. 

5.2.1 Electric Engines 

Electric engines provide clean power. Since no combustion takes place in these 
systems, noise pollution is minimized. Furthermore, electric engines are easier to 
maintain when compared to internal combustion engines. Starting an electric engine 
simply requires a DC power source rather than the fuel, glow plugs and starters needed 
for internal combustion engines. Hence, electric engines need less time and effort for 
starting, steady operation and maintenance. 
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However, the electric engine does have some considerable disadvantages. The 
most significant is its low thrust to weight ratio (T/W). Batteries, the clean power source 
for the electric engine, are also its primary detriment. For example, an electric engine 
which weighs 4.39 N requires three standard 6 Volt batteries weighing 8.01 N. This 
engine’s maximum thrust output is 15.57 N, yielding a max T/W of only 1.26. A 
comparable internal combustion engine, such as the one chosen for the flying wing, has 
a max T/W using manufacturer’s data of 2.43. A higher T/W allows the flying wing to 
have better maneuverability in the air, and ensures enough reserve of thrust to maintain 
flight. 

5.2.2 Internal Combustion Engines 

In the discussion of the electric engines many of the disadvantages of the internal 
combustion engines were mentioned. These engines produce plenty of pollution. The 
testing and operation of these engines must be done in an environment where fresh air 
can be circulated. Also, two-cycle engines produce tremendous noise pollution, although 
a muffler attached to the engine will reduce some of the noise pollution. Furthermore, 
air pollution is created by the exhaust. The unburned castor oil (required for lubrication) 
spills out from the exhaust, creating further pollution. 

The internal combustion engines must be maintained to ensure smooth operation 
and reduce wear and tear due to the combustion process. Great care is required to assure 
that no foreign particles enter the combustion chamber. The engine must be cleaned 
externally, so that debris do not collect on residual oil that lingers on the engine’s body. 
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Some hands on experience is required in the engine’s operation. Choosing an 
optimum fuel/air mixture and controlling the throttle during testing requires careful 
percision. 

The main advantage of the internal 
combustion engine is its high thrust to 
weight ratio, since fuel provides more energy 
per density compared to available batteries. 

For this reason, the internal combustion 
engine was chosen for the flying wing. 

Furthermore, internal combustion engines 
are widely accepted as the norm in the model 
airplane industry with the exception of sail 
planes. 

The flying wing is powered by two 4.07 cc O.S. VF-DF Engines (Figure 5.1). 
Each engine produces 820.30 Watts at 22,000 rpm and weighs 2.23 N. 

5.3 Method of Propulsion 

Two main methods of propulsion were considered, propellers and ducted fans. 
The first method is widely used and accepted in the model airplane industry. The second 
is a rather new approach and less accepted mode of propulsion. 



Figure 5.1: 4.07 cc O.S. Engine 
(O.S. Engines Catalog . 1990) 
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5.3.1 Propellers 

Free propeller driven airplanes are the norm in the model airplane industry. They 
have been field tested quite extensively and have proven to be reliable. The propulsion 
system consists of a propeller, an engine and the fuel tank. The propeller can be easily 
mounted to the driving shaft of the engine. There are no requirements for other major 
accessories, resulting in reduced weight. Furthermore, the engine is readily accessible, 
since it is not encased in a duct, in contrast to the ducted fan setup. 

There are two major disadvantage with a free propeller setup for the flying wing: 
propeller wash which interferes with control surfaces, and the large ground clearance 
required for a large diameter propeller. 



all control surfaces clear of any disturbances, such as those created by propeller wash. 

Since propellers have a relatively large diameter compared to the thickness of the 
flying wing, large landing gear would be required to attain the needed ground clearance. 
Utilizing large landing gear adds unnecessary drag. However, minimumizing drag was 
a design goals of the aerodynamics team. Therefore, to maintain stability and 
aerodynamic integrity, a conclusion was reached not to use a propeller. 
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5.3.2 Ducted Fan Engines 

Although ducted fan engines have a much more complex assembly process than 
regular propeller engines, their clean aerodynamic shape seemed to fit perfectly in the 
airfoil, producing enough thrust and yet keeping drag to a minimum. 

The ducted fan engine consists 
of a piston engine enclosed in a 
cylindrical, duct, with a small seven- 
blade rotor spinning at much higher 
rpm’s than regular propeller engines to 
produce the same amount of thrust 
(Figure 5.3). Since the engine is 
completely enclosed inside the 
cylindrical housing, it is very hard to 
assemble the kit or to do any 
maintenance on the engines. In fact if any problems arise with the engines, the entire 
kit would have to be removed from the aircraft and then disassembled. The assembly 
process is rather complicated because extremely tight tolerances have to be kept in order 
for the fan to work properly. The rotor has to be perfectly balanced, since it rotates at 
speeds around 22000 rpm. The slightest imbalance would generate undesired vibrations. 



Figure 5.3: Ducted Fan Engine Layout 

(RK:.7. 20 MK III Instructions and 

Specifications . 1992, page 3) 
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5.4 Propulsion Systems Configuration 



Figure 5.4: Isometric View of Propulsion System Layout 


The final configuration of the propulsion system is shown in Figure 5.4. The 
three dimensional drawing clearly shows the shell, the rotor, the engine, the stator and 
all the supports. The engine’s carburetor and cylinder head protrude from the duct. Each 
engine assembly as shown, weighs 4.00 N. Both engines receive fuel from a central fuel 
tank located near the center of gravity of the airplane. The manufacturer’s ducted fan 
test results are shown in Figure 5.5. The operating speed is approximately 22000 RPM, 
with a resulting thrust of 15.60 N for each engine, which provides more than enough 
combined thrust to propel the flying wing. 

The experimental thrust obtained from strain-gauge measurements is lower than 
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Figure 5.5: Thrust vs. RPM 

the manufacturer’s results. At 22000 RPM, only 8.9 N of thrust is acquired (Figure 
5.5). Even with this lower value, the combined two engine thrust of 17.8 N is more 
than sufficient to overcome the empirical drag calculated by the aerodynamics sub-group. 
The ultimate test of the propulsion system will come during the ground and flight testing 
of the flying wing. 
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CHAPTER SIX: CONCLUSIONS AND 
RECOMMENDATIONS 


Many engineering difficulties were encountered during the design and construction 
of the flying wing, Elang. After the plane construction was completed, there appeared 
to be many components and processes which could be further optimized through more 
research, development and testing. Of course many of these revelations were not 
obvious to the project team before the actual construction began. The performance of 
Elang depends on the following criteria: overall efficiency of the propulsion system, 
structural design, material selection, stability, aerodynamic analysis and the overall 
weight of the plane. 

Even though, the propulsion system generated an adequate amount of power, it 
seemed to generate quite a bit of problems. Since this was original thought as an oblique 
flying wing two engines were used. However the ability to synchronize the power 
produced by the engines became more difficult then anticipated. Though the engines 
preformed fairly well at higher revolutions per minute, at lower revolutions per minute 
the performance declined. The engines would simply stall. Another point of difficulty 
was the length of the fuel lines. Each engines draws fuel over a 25 cm line. This 
appeared to be too much for the small 0.25cc engines to handle. The inability to place 
a muffler or a tune pipe posed a problem for the design team. For a future design it is 
recommended to use one large engine or to design an adequate throttle control. A good 
consideration for throttle control could be one that measures the revolutions per minute 
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in flight and adjusts for any deviations. The incorporation of tune pipes or mufflers 
could prove extremely pleasant to the ears. 

Stability was of course one of the primary problems with a plane of this type. 
Finding the center of gravity and designing around that proved extremely difficult. Also, 
the sizing of the multiple control surfaces was extremely important. The group decided 
to oversize all control surfaces to account for the inherent instability of an aircraft 
without a fuselage or horizontal stabilizer. Finally, an airfoil was designed with a 
positive moment about the aerodynamic center. This characteristic is essential in 
counteracting the inherent instability of the flying wing. 

Aerodynamically, the project group found that the use of computer codes for 
analysis, although fast and relatively easy, proved to be quite risky. It is recommended 
that the computer codes used should be verified for proper accuracy and operation prior 
to using them for design and analysis. Also, not being able to build a prototype model 
for wind tunnel testing limited the group to purely theoretical aerodynamic analysis. 

The aircraft uses three landing gears. One steerable nose gear and two back 
landing gears. Even though, the design of the landing gears was left for the last minute 
they seem to preform quit well. The only minor problem was the length and width of 
the axle which held the front gear. However once the length was shortened and the 
diameter of the axle increased the problem was solved. 

Structurally, Elang was very solidly built. The carbon fiber and vacuum bagging 
technique proved to be very effective. The foam core, although difficult to shape 
correctly the first time, became quite easy to work with and proved to be the best mold 
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for the carbon fiber skin. However, the number of carbon composite layers used around 
the engine, could be reduced from four to maybe two or three layers for simpler 
construction. 

The participation of outside vendors in this project assisted the building team. 
The use of a 2-D CAD driven hot wire at Foam Technology, saved valuable time and 
minimized the amount of sanding needed by hand. This produced a more accurate mold 
which in turn provided a better surface for the carbon to lie on. 

Finally, in the area of group dynamics, this project group found that thirteen 
group members seemed to be too large. It is recommended that, if at all possible, a 
smaller group size of approximately eight people would be best. This way, 

communication and coordination of activities will be that much easier. 

Overall, the project team was very satisfied with the analysis, design, 
construction, and performance of Elang. The recommendations mentioned above indicate 
areas in which the project team felt limited. Most of these recommendations were 
realized through experience and did not come to light until the completion of the aircraft. 
Further research and development in these areas are encouraged because the possibilities 
for various design configurations of this type of aircraft are numerous. 
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APPENDIX A 


MATHCAD RESULTS FOR STABILITY ANALYSIS 
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Read Pints from data file (hp33.pm-section points, hp33int.pm-lift, drag, moment) 


i =0-45 j =0.. 12 

point - RE ADPRN( points) LDM = READPRN(Idm) 

x a i := P° im i + 45,o yu. ;= P° int i+ 45.i >1. r P° int 45-i,i 
a 0 : = LDM j 0 C j := LDM. j C d .:=LDM j>2 C 4 
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Spline fit th e lift, drag and moment curves 
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Design Parameters 

Weight wing 

Weight engine empty 
Weight engine with fuel 

Wing x center of gravity 
- 0.0056 l-m Wing y center of gravity 

Cruise velocity 

Cruise thrust 


w w 16-N 


Wee :=5.719-N 


w ef =9.578-N 


x wcg ;=0 - 209 rT1 


m 
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T e :*».N 
e 30 
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2 
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y q) := 0 -m qs=q-S w e :=w 


Dynamic pressure 
Wing surface area 


ee 


P : = 0 -deg 
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p := 1 . 225 -— 
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Sweep Angle 
Engine pitch angle 

Sea level density 


«n :=4 
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C L ( a n) =0-536 
C D ( a n) = 8.91 i* 10 

C m4 (a n ) =0011 


“4 


AR -1 


Spline airfo il coordinates 

vY u := cspline^x a -cord, y u -cordj 
vY | :=cspline^x a -cord,y j-cordj 


Y u ( x) : = interp^vY u , x a -cord, y u -cord , xj 

Y j(x) :=interp(vY j.Xa-cord.y j-cord.xj 


Aerodynam ic center, and C ]yf ac 


qr 


^ mac^ ‘ 


/ 1 C m4 ( a)\ 

i- -cord x. c (a) := 

l-Ctfa) 

X (a) + C L (a)-— X_(a) 

\ 4 C L (a) j 

Ida j 

da 


cord 


da 


( x ac( a ) ^ cp^^) 


a:=- 2 f - 1 . 5..10 



0.124 


0.025 


^mac vs - a 


X ac vs. a 





0.022 

0.019 

0.016 

0.013 

0.01 
















v 





\ 








~2 1.2 4.4 7.6 10.8 14 

a 


^ mac^ C mac( a )- cos (P> 
M ac (a ^ :=C mac (a) -<l s - cord 


"2 1.2 4.4 7.6 10.8 14 

a 


mac 




0.019 


M ac ( a r) = 1 088* m *N 


C Lu( a ) : =C L (a)-cos(a-deg) 
C Lb ( a ) : = C L (a)-sin(a-deg)-cos(( 3 ) 

L u (a) : = C j^/oO-qs 
L b (a) =C Lb (a) qs 

L u( a J = 65.769 *N 
L b (a n ) = 4.599 *N 

:=x ac (a) 

T e20 =T e' cos( P) 

T f : " T e20 cos(r) 

T f = 1.667 -N 


C Db (a) =C D (a)-cos(a-deg)-cos(P) 
c Dd^ :=c D(tt)' sin (° t 'd e g) 

D b (a) = C Db (a)qs 
D d (a) = C Dd (a)qs 

D d( a n) =-0-008 *N 
D b( a n) =0-109*N 


T e20 = 1.667-N 

T d :=T e20' sin(r) 
T d = 0 *N 


Extra mass estimation 


w mass ' 0.5-kg-g 

x mass ’ y mass - y WC g 

x (\ W* ecg ' 

(^‘ w e) + x wcg' w w + x mass‘ w mass 

cgl, x ecgj ' 

( 2 ' w e + w w + w mass) 

v (v V- Yecg ' 

(2-w + y wcg’ w w f y mass' w mass 

ycgi y ecgj ■ 

( 2 ' w e+ w w+ w mass) 

^( x ecg ,y ecg ,<x ) 

= M ac (a) ... 

+ L u ( ct)- (x C g( x eC g) _ x l(®)^ ••• 

+ -L b ( a )-(y C g(yecg) ~ y cp) ■" 

+ D d ( a )-(x C g(x eC g} - x L (a)) ... 

+ -D b (a)-^y C g^y ecg) - yep) 




^ M( x ecg ,y ecg ,a ) ^'mac^ a ^'" 

+ C L u (a)-(x cg (x ecg ) - x L^ a) )'^^ - 

+ - C Lb( “) • (y cg (y ecg ) - Y C p) *“ - 

+ C Dd (a)-(x cg (x ecg )-x L (a))--^... 

+ -C Db( a )'(y C g(Yecg) _y cp^’^ 

Solve The above equations in order to detremine the position of the engines 

a n =4 Y ecg := y wcg x ecg :=0 ’ m x ecg :=root ( M ( x ecg’ y ecg’ a n) ,x ecg) 

x ; = -0.1*m,0*m.. cord 





Moment About the Center of Gravity Solution: 

w e = 5.719’N 

x ecg =-0.01657*m x cg (x ecg ) = 0.10512-m 

yecg =0.00561-m y cg ( yecg ) -0.00561-m 

w mass =4 - 9 °3*N 

-o.l 0 0.1 0.2 0.3 0.4 

X x mass = 005 * m y mass = 0 006 ' m 

a :=-2,- 1.5.. 10 

Static Stability; Coefficent of Moment About the Center of Gravity 

0.027 
0.016 

r ( „ 0.006 
C M \ x ecg ■ y ecg * °y 

0 -0.005 

“0.015 
-0.026 

-2 0 2 4 6 8 10 




a 





x-hb yc(x) -O.OOx 


xx - 0*m, 0.0 1-m.. cord 


.02*m x f ~ x ecg f 0.02-m x = x j ,x j + 0.01*m.. x f 

Side View 



- 1.25 m,- 1.20-m.. t.25m 


y e i< x) 



Tod View 



x.x.x.O-m.xx.xx 








vC Ldef =cspline(a D0 ,C Liefo) C Ldef (a) := inter P( vC Ldef’ a DO' C LdefO’ a ) 

vC Ddcf : = cspline(a D0 ,C Ddef0 ) C D def (a) " imer P( vC Ddef’ a D0> C DdefD’ a ) 

vC C4def := cs P line ( a D0 , C Ddcf0 ) C C 4def (a) =inter p(vC C4def’ a D0' C C4def0’ a ) 

a q - 30.. 30 

vs. a Cq vs. Cl 




c Ldef( a d) 


^m4 vs - a 

1000 
400 

C Ldef( a p) -200 
c Ddef( a d) -800 
-1400 
-2000 

-30 “18 -6 6 18 30 

a D 


Cl/Cq vs. a 





i del 

s^,f =0.02S s Hpf =0.018*m width = width =0.21 8*m 

aet .18-cord 


qs def r Q' s def qs def = 2.461 ‘N qs w :=qs - qs def qs w = 1 20.605 *N 

C Ludef( a d) : = C Ldef( a d) ' cos ( a D' de g) L udef( a d) " C Ludef( a d) ’V def 

C Lbdef( a d) ; = C Ldef( a d) - sin ( a D' de s) cos( P> L bdef( a D) : " C Lbdef( a d) def 

C Dbdef( a d) := C Ddef ( a d) - cos ( a D‘ de g) - cos < P) D bdef( a d) = C Dbdef( a d) ' V def 

c Dddef( a D) :=C Ddef( a d) ' sin ( a D' de g) D ddef(« d) :=C Dddef( a d) ^ def 


M C4def( a d) >C C4def( a d) '^ s de f cord-cos ( P) 






control Surfaces 

Read Pints from data file (flap33.pm-lift, drag, moment) 

j =0.. 32 

LDM d : = READPRN( flap) 


a DO. :_LDM d 


j.o 


C LdefU ~ LDM d 

J j.> 


C DdefO. LDM d 


j.2 


C C4defO. :=LDM d. , 

J J.6 


Angle of Attack 


Lift Coefficient Drag Coefficient Moment Coefficient 

about quater chord 


-30 


2.295 ' 


0.00664 


-5.021 

-26 


2.074 


0.00624 


-3.804 

-22 


1.847 


0.00557 


-2.757 

-18 


1.615 


0.00472 


-1.857 

-16 


1.498 


0.00426 


-1.464 

-12 

C LdefO “ 

1.261 

C DdeFO = 

0.0033 

C C4defO = 

-0.801 

-10 


1.141 


0.00283 


-0.552 

-9 


1.081 


0.0026 


-0.446 

-8 


1.021 


0.00238 


-0.35 

-7 


0.961 


0.00216 


-0.264 

-6 


0.9 


0.00195 


-0.19 

-5 


0.84 


0.00175 


-0.133 

-4 


0.779 


0.00155 


-0.085 

-3 


0.718 


0.00137 


-0.052 

-2 


0.658 


0.0012 


-0.03 

-1 


0.597 


0.00104 


-0.015 

0 


0.536 


0.00089 


0.011 

1 


0.475 


0.00076 


0.05 

2 


0.414 


0.00063 


0.103 

3 


0.353 


0.00052 


0.167 

4 


0.292 


0.00043 


0.245 

5 


0.231 


0.00035 


0.334 

6 


0.169 


0.00028 

* 

0.437 

7 


0.109 


0.00023 


0.55 

8 


0.048 


0.00019 


0.676 

9 


-0.013 


0.00017 


0.814 

10 


-0.074 


0.00016 


0.965 

12 


-0.196 


0.00017 


1.302 

16 


-0.437 


0.00035 


2.107 

18 


-0.557 


0.0005 


2.571 

22 

j 

-0.794 


0.00085 


3.611 

26 


-1.027 


0.00118 


4.783 

30 . 


-1.256 


0.00138 


6.068 



li i 1 1 ubliiil 


a Dn =4 

L u( a d) :=C Lu( a d) -c 1 s w L b( a d) : = C Lb( a D)‘^ s w 
D d( a d) ’ C Dd( a d) '^ s w D b( a d) = C Db( a d) w 

M ac ( a d) : = C mac( a d) W w' cos ( P> cord 

^ def( x ecg’ y ecg’ ® D ~^ac^^ ■” 

+ L u ^ a ^'( x cg( x ecg) “ x l) 

+ - L b( °t) • (y C g(y ecg) _ y cp) ■" 

+ D d^a)'( x C g( x ecg) _ x l) — 

+ - D b (a)-(y c „(y ecg ) - y cp ) - 

+ M C4def( a Dj - 

+ ^ udef( a d) ' ( x cg( x ecg) _ x Ldef) — 
+ ’ ^ bdef ( a d) ' (y cg(y ecg) ~ y cp) ••• 
+ ^ ddef( a d) '( x cg( x ecg) ~ x Ldef) •" 
+ " ® bdef ( a d) ' (y cg(y ecg ) “ y cp) 


Moment About the Center of Gravity (For Control Surface) 





Angular Acceleration 


Angular Acceleration 


angular a a ( a D‘ a ri) 


200 

120 

40 


rad 


“120 

- 200 ' 






Lz 





/ 




T 7 ^ 







/ 

r 




7^ 






"30 “18 “6 6 18 30 

a n 



Aileron Analysis 


s def :=0.02-S s def = 0.018*m 2 width 

S Hef 

;= u width -0.21 8* m 

0.1 8 -cord 

9 s def : = c l' s def ^ w : = qs - qs def 

qs w = 120.605 *N 

C Ludef (« d) : = C Ldef(“ d) - cos ( a D' de g) 

^udef( a d)- c Ludef ( a d) ‘^ s def 

C Dddef( a d) ~~ C Ddef ( a d) ‘ sm ( a D' de §) 

D ddef( a d) = C Dddef( a d) V s def 

L u(“d) :rC Lu( a D)-q s w 


^ ldef ( a d) - ^ udef( a d) 
^ rdef ( a d) ' ” ^ udef( a d) 

D ldef ( a d) ;=D ddef( a d) 
^ rdef( a d) ^ ddef( a d) 

x ccg "°- m y ecg :=0 ' m 

I ^1= 0.62161 1-kg-m 2 

z al '~ ^ ' m z ar~ z al 


M rdef( a d) : = L ldef ( a D + 30 ) z al - 


+ ' rdef{ 3 ® “ a d) z ar) 

+ ' ( D ldef( a D + 30 )' z al) - 

+ D rdef{ 30 _ a d) ' z ar 


, x M rdef(“ d) 
angidar a a( a Dy := 

1 XX 



XX 




Moment About the Center of Gravity (For Control Surface) 


Angular Acceleration 


M rdcf 




Rudder Analysis 

Read Pints from data file (sd 8020 .pm-lift, drag, moment) 
j := 16.. 32 

LDM r : = READPRN( rudder) 


a 


rl. 


:= LDM , 


a.i = 


j- 16 

Angle of Attack 
0 
1 
2 

3 

4 

5 

6 

7 

8 

9 

10 
12 
16 
18 
22 
26 
30 

h : = 0.. 16 


j.o 


•Lr. 


:= LDM 


d. 


a r0 a rl 
h 


C T . = 


j— 16 j — 16,1 

Lift Coefficient 

2.295 
2.074 
1.847 
1.615 
1.498 
1.261 
1.141 
1.081 
1.021 
0.961 
0.9 
0.84 
0.779 
0.718 
0.658 
0.597 
0.536 


Dr. 


: = LDM , 


DrO : ~ c 

n 


a r0 '~ a rl 

( 16— h) ( 16-j-h) <h) 

Dr,, ^ DrO.,, 

(16-h) (16-t-h) 


c Dr = 


•j-16 j,2 

Drag Coefficient 

0.00089 
0.00076 
0.00063 
0.00052 
0.00043 
0.00035 
0.00028 
0.00023 
0.00019 
0.00017 
0.00016 
0.00017 
0.00035 
0.0005 
0.00085 
0.00118 
0.00138 


LiO. 


:=-c , 


’ C4r. 


LDM 


d. 


j-16 j > 6 

Moment Coefficient 
about quater chord 

0.011 
0.05 
0.103 
0.167 
0.245 
c C4r = | 0.334 
0.437 
0.55 
0.676 
0.814 
0.965 
1.302 
2.107 
2.571 
3.611 
4.783 
6.068 


‘(16-h) 


c n :=Cc4r (i6-h> ( 


C4r0 


(16+h) 


Dr, 


LrO 


•-Ci 


{ 16-j-h) 





Spline fit the lift, drag and moment curves for rudder control 


vC Lr -cspline^ro.CLrt)) 

a) : = interp^vC Lr Q.ct) 

vC Dr ; =csplineia l0 .C Dr0 J 

C o r (a) -interpjvC Or’ a r0’ C Dr0’ a ) 

C4r -csplinc^a^.C p,^) 

C C4r ( a) := inteip(vC C4def , a *), C C4r0 , a) 

cord r = 0.15-m 



S r :=0.02S qs r : = q-S r 

Lenght r : = 0.3 -cord r 

S r 

width r : = S rc width . 

r cord r rs 

S r = 0.018*m 2 qs r = 2.461 *N 

Lenght r =0.04 5* m 

width r =0.119*m S rs =0.005* 

L r( a r) :=c Lr( a r) -< l s r L rz 

(a r ) :=L r (a r )-cos(p) 

L rx(“r) : "L r (a r )-sin(p) 

D r(“r) : =CDr(a r )*qs r D n 

'( a r) : =D r (a r )-sin(P> 

D rx( a r) := D r (a r )-cos(P) 

cord r 



x Lr •= 4 x Lrz =: 

x-^-sin(P) XLn4 : 

= x Lr'COs( P) 


Z| r :=0.57-m X| r -0.1-m 

z^^O^-m x rr :=0.1-m 


M lr( a r) ( L rx( a r) + D rx( a r)) ‘( z lr x Lrz) ( L rz( a r) + ^ rz( a r)) '( x lr x Lrx) 

M n( a r) ' - " [ rx ( a r) + ^ rx ( a r)) ' ( z rr + x Lrz) ] ~ rz( a r) D rz( a r)) ' ( x rr +" x Lrx) 

M r( a r) :=M lr(«r) ^ M rr( a r) 

a r :=- 30.. 30 


Moment About the Center of Gravity (For Control Surface) 

width r = 0.119 th 
L enght r = 0.045* m 
s rs = 0.005* m 2 

-30 -18 -6 6 18 30 

«r 



I yy :=0.656673-kg-m 2 


angular ^o,) 




yy 





APPENDIX B 


COORDINATES AND DATA FOR HPRS 33 
AIRFOIL SECTION 


hi 


Coordinates for HPRS 33 Airfoil Section 


X 

0.00000 
0.00182 
0.00455 
0.01136 
0.02273 
0.03409 
0.04545 
0.06818 
0.09091 
0.11364 
0.13636 
0.15909 
0.18182 
0.22727 
0.27273 
0.31818 
0.36364 
0.42209 
0.46755 
0.51300 
0.53573 
0.55845 
0.58118 
0.60391 
0.62664 
0.64936 
0.67209 
0.69482 
0.71755 
0.74027 
0.76300 
0.78573 
0.80845 
0.83118 
0.85391 
0.87664 
0.89936 
0.91559 
0.93182 
0.94319 
0.95455 
0.96591 
0.97727 
0.98864 

1.00000 

1.00000 

l. 00000 


Y(upper) 

0.00000 

0.01182 

0.01850 

0.02790 

0.03786 

0.04522 

0.05091 

0.05965 

0.06645 

0.07190 

0.07648 

0.08044 

0.08372 

0.08889 

0.09245 

0.09463 

0.09545 

0.09505 

0.09335 

0.09015 

0.08770 

0.08522 

0.08194 

0.07822 

0.07404 

0.06945 

0.06451 

0.05925 

0.05375 

0.04745 

0.04222 

0.03621 

0.03012 

0.02399 

0.01875 

0.01363 

0.00951 

0.00700 

0.00500 

0.00425 

0.00375 

0.00330 

0.00415 

0.00550 

0.00700 

0.00500 

0.00300 


Ydower) 

0.00000 

-0.00860 

-0.01275, 

-0.01812 

-0.02377 

-0.02818 

-0.03152 

-0.03717 

-0.04151 

-0.04493 

-0.04792 

-0.05034 

-0.05134 

-0.05234 

-0.05234 

-0.05234 

-0.05234 

-0.05234 

-0.05234 

-0.05234 

-0.05134 

-0.05019 

-0.04766 

-0.04470 

-0.04132 

-0.03767 

-0.03383 

-0.02988 

-0.02588 

-0.02174 

-0.01793 

-0.01402 

-0.01052 

-0.00752 

-0.00550 

-0.00466 

-0.00466 

-0.00475 

-0.00420 

-0.00400 

-0.00380 

-0.00300 

- 0.00100 

0.00100 

0.00700 

0.00500 

0.00300 



APPENDIX C 


COORDINATES FOR SD 8020 AIRFOIL SECTION 


IV 



Coordinates for SD 8020 Airfoil Section 


X 

Y 

*■ upper 

0.0000 

0.0000 

.00276 

.00645 

.01065 

.01345 

.02318 

.02041 

.04024 

.02697 

.06179 

.03287 

.08774 

.03802 

.11789 

.04233 

.15203 

.04574 

.18987 

.04824 

.23106 

.04982 

.27523 

.05051 

.32196 

.05034 

.37077 

.04938 

.42116 

.04770 

.47261 

.04536 

.52456 

.04246 

.57647 

.03909 

.62777 

.03535 

.67789 

.03135 

.72627 

.02722 

.77237 

.02308 

.81560 

.01908 

.85538 

.01532 

.89118 

.01188 

.92247 

.00876 

.94885 

.00591 

.97017 

.00330 

.98625 

.00131 

.99646 

.00027 

1.00000 

0.00000 


X 

Y 

A lower 

0.0000 

0.00000 

.00276 

-.00645 

.01066 

-.01345 

.02319 

-.02041 

.04024 

-.02697 

.06180 

-.03287 

.08774 

-.03802 

.11790 

-.04233 

.15204 

-.04574 

.18987 

-.04824 

.23107 

-.04982 

.27524 

-.05051 

.32197 

-.05034 

.37077 

-.04938 

.42117 

-.04769 

.47262 

-.04536 

.52457 

-.04246 

.57648 

-.03908 

.62778 

-.03534 

.67790 

-.03135 

.72629 

-.02722 

.77238 

-.02308 

.81561 

-.01907 

.85539 

-.01532 

.89119 

-.01188 

.92248 

-.00876 

.94886 

-.00591 

.97018 

-.00330 

.98626 

-.00131 

.99647 

-.00027 

1.00000 

0.00000 



n 


APPENDIX D 


AERODYNAMIC DATA FOR AIRCRAFT 


v 


w»a a i r 

densi ty 

Temp Cruis 

R 

Speed of Sound 


< 1 bm/f t*3> 

<F) 

<f t*lbf/lbm«R> 

(ft/s) 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

l .4 

0 .0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0 .0764743 

59.02 

53.34 

1 116.792264 

1 .4 

0 .0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0 .0764743 

59.02 

53.34 

1116.792264 

l .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1.4 

0 .0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 116.792264 

1 .4 

0 .0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1116.792264 

1 .4 

0 .0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0 .0764743 

59.02 

53.34 

1 1 16.792264 

l .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0 .0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1116.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1116.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 116.792264 

1 .4 

0 .0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0 .0764743 

59.02 

53.34 

1116.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0 .0764743 

59.02 

53.34 

1116.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1116.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1.4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1116.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1116.792264 

1 .4 

0.0764743 

59.02 

53.34 

1116.792264 

1 .4 

0.0764743 

59.02 

53.34 

1116.792264 

1.4 

0.0764743 

59.02 

53.34 

1116.792264 

1.4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 

1 .4 

0.0764743 

59.02 

53.34 

1 1 16.792264 



V Cru i se 

V Cruise 

Mach Cruise 

Beta 

Cl aloha 


(m/s) 

<• ft/s) 



sect i on 


0 

0 

0 

1 

6.97 


4 

13.12335958 

0.011750940621 

0.99986191539 

6.97 


8 

26.24671916 

0.023501881241 

0 .99944766158 

6.97 


12 

39.37007874 

0.035252821862 

0.99875723855 

6.97 


16 

52.49343832 

0.047003762483 

0.99779064631 

6.97 


20 

65.6167979 

0.058754703103 

0 .99654788486 

6.97 


0 

0 

0 

1 

6.97 


4 

13.12335958 

0.011750940621 

0.99986191539 

6.97 


8 

26.24671916 

0.023501881241 

0.99944766158 

6.97 

— 

12 

39.37007874 

0.035252821862 

0.99875723855 

6.97 


16 

52.49343832 

0.047003762483 

0.99779064631 

6.97 

- ■ 

20 

65.6167979 

0.058754703103 

0.99654788486 

6.97 

^ , 

0 

0 

0 

1 

6.97 


1 

3.280839895 

0.002937735155 

0.99999136971 

6.97 

- . 

2 

6.56167979 

0.00587547031 

0 .99996547885 

6.97 


3 

9.842519685 

0.008813205465 

0.99992232741 

6.97 


4 

13.12335958 

0.011750940621 

0.99986191539 

6.97 


5 

16.404199475 

0 .014688675776 

0.9997842428 

6.97 

■ — * 

6 

19.68503937 

0.017626410931 

0.99968930964 

6.97 

E3 

7 

22.965879265 

0 .020564146086 

0.9995771159 

6.97 


8 

26.24671916 

0.023501881241 

0.99944766158 

6.97 

* — 

9 

29.527559055 

0 .026439616396 

0.99930094668 

6.97 


10 

32.80839895 

0.029377351552 

0.99913697122 

6.97 


1 1 

36.089238845 

0 .032315086707 

0.99895573517 

6.97 


12 

39.37007874 

0.035252821862 

0.99875723855 

6.97 


13 

42.650918635 

0.038190557017 

0.99854148135 

6.97 

w 

14 

45.93175853 

0.041128292172 

0.99830846358 

6.97 


15 

49.212598425 

0 .044066027327 

0.99805818524 

6.97 

- ■ 

16 

52.49343832 

0.047003762483 

0.99779064631 

6.97 


17 

55.774278215 

0.049941497638 

0.99750584681 

6.97 


18 

59.0551181 1 

0.052879232793 

0.99720378674 

6.97 

“ 

1? 

62.335958005 

0 .055816967948 

0.99688446609 

6.97 


20 

65.6167979 

0.058754703103 

0.99654788486 

6.97 


0 

0 

0 

1 

6.97 


4 

13.12335958 

0.011750940621 

0.99986191539 

6.97 

= 

8 

26.24671916 

0.023501881241 

0.99944766158 

6.97 


12 

39.37007874 

0.035252821862 

0.99875723855 

6.97 


16 

52.49343832 

0 .047003762483 

0 .99779064631 

6.97 


20 

65.6167979 

0 .058754703103 

0.99654788486 

6.97 


0 

0 

0 

1 

6.97 


4 

13.12335958 

0 .01 1750940621 

0.99986191539 

6.97 

- 

8 

26.24671916 

0.023501881241 

0.99944766158 

6.97 

__ 

12 

39.37007874 

0.035252821862 

0.99875723855 

6.97 


16 

52.49343832 

0.047003762483 

0.99779064631 

6.97 

_. - 

20 

65.6167979 

0.058754703103 

0.99654788486 

6.97 


4 

13.12335958 

0.011750940621 

0 .99986191539 

6.97 

*- 

8 

26.24671916 

0.023501881241 

0.99944766158 

6.97 


12 

39.37007874 

0.035252821862 

0.99875723855 

6.97 


16 

52.49343832 

0.047003762483 

0.99779064631 

6.97 


20 

65.6167979 

0.058754703103 

0.99654788486 

6.97 



Sref 

C root 

Sexposed Aspect 

effSweep t/c 

CL aloha 

Alpha not 

<f t*2) 

(ft) 

<f t*2> 


max 


( deg) 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.071 1746426 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0716825682 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0732069238 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0757494475 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0793130414 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0839017787 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.071 1746426 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0716825682 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0732069238 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0757494475 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0793130414 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0839017787 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.071 1746426 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0712063851 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0713016149 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0714603388 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0716825682 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0719683187 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0723176108 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0727304693 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0732069238 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.073747008 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0743507606 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0750182246 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0757494475 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0765444817 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0774033837 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.078326215 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0793130414 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0803639335 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0814789665 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0826582201 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0839017787 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.071 1746426 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0716825682 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0732069238 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0757494475 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0793130414 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0839017787 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.071 1746426 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0716825682 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0732069238 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0757494475 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0793130414 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0839017787 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0716825682 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0732069238 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0757494475 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0793130414 

1 

9.6106343 

1 .46 

9.6106343 

7 

0 

5.0839017787 

1 


Cl not 

Cru i se 

Alpha crui 

** CL »* 

q»S 

L 

L 


ancjl e 

< rad > 

< wi ng > 

( 1 b-f ) 

wing 

< N > 

0.081 

0 

0 

0.081 

0 

0 

0 

0.081 

0 

0 

0.081 

1 .9654941 

0 .159205 

0.7081792 

0.081 

0 

0 

0.081 

7.8619765 

0.6368201 

2.8327169 

0.081 

0 

0 

0.081 

17.689447 

1 .4328452 

6.3736131 

0.081 

0 

0 

0.081 

31 .447906 

2.5472804 

1 1 .330868 

0.081 

0 

0 

0.081 

49.137353 

3.9801256 

17.704481 

0.081 

2 

0 . 0349066 

0.25801738891 

0 

0 

0 

0.081 

2 

0 .0349066 

0.25803511886 

1 .9654941 

0 .5071665 

2.255989 

0.081 

2 

0 .0349066 

0.25808832891 

7.8619765 

2.0290844 

9.025817 

0.081 

2 

0.0349066 

0 .25817707973 

17.689447 

4.5670098 

20.315072 

0.081 

2 

0 .0349066 

0.25830147262 

31 .447906 

8.1230405 

36.133084 

0.081 

2 

0 .0349066 

0 .25846164977 

49.137353 

12.700121 

56.492955 

0.081 

4 

0.0698132 

0.43503477783 

0 

0 

0 

0.081 

4 

0.0698132 

0 .43503699387 

0.1228434 

0.0534414 

0.2377193 

0.081 

4 

0 .0698132 

0.43504364217 

0.4913735 

0.2137689 

0.9508916 

0.081 

4 

0.0698132 

0.43505472319 

1.1055905 

0.4809923 

2. 1395606 

0.081 

4 

0.0698132 

0.43507023772 

1 .9654941 

0.855128 

3.8037989 

0.081 

4 

0.0698132 

0.43509018687 

3.0710846 

1 .3361988 

5.9437082 

0.081 

4 

0.0698132 

0.43511457206 

4.4223618 

1 .9242341 

8.5594196 

0.081 

4 

0.0698132 

0.43514339503 

6.0193258 

2.6192699 

1 1 .651093 

0.081 

4 

0.0698132 

0 .43517665782 

7.8619765 

3.4213487 

15.218917 

0.081 

4 

0.0698132 

0.43521436281 

9.9503141 

4.3305196 

19.2631 1 1 

0.081 

4 

0.0698132 

0.4352565127 

12.284338 

5.3468383 

23.783922 

0.081 

4 

0.0698132 

0.43530311047 

14.864049 

6.4703669 

28.781626 

0.081 

4 

0.0698132 

0.43535415946 

17.689447 

7.7011744 

34.256531 

0.081 

4 

0.0698132 

0 .43540966332 

20.760532 

9.0393362 

40.208971 

0.081 

4 

0.0698132 

0 .43546962599 

24.077303 

10.484934 

46.63931 1 

0.081 

4 

0.0698132 

0.43553405176 

27.639761 

12.038057 

53.547946 

0.081 

4 

0.0698132 

0.43560294525 

31 .447906 

13.698801 

60.935301 

0.081 

4 

0.0698132 

0 .43567631 136 

35.501738 

15.467266 

68.801828 

0.081 

4 

0.0698132 

0.43575415534 

39.801256 

17.343563 

77.148011 

0.081 

4 

0.0698132 

0.43583648278 

44.346461 

19.327806 

85.974363 

0.081 

4 

0.0698132 

0.43592329955 

49.137353 

21 .420117 

95.281428 

0.081 

6 

0 .1047198 

0 .61205216674 

0 

0 

0 

0.081 

6 

0 .1047198 

0 .61210535658 

1 .9654941 

1 .2030895 

5.3516087 

0.081 

6 

0 .1047198 

0 .61226498673 

7.8619765 

4.813613 

21 .412017 

0.081 

6 

0 .1047198 

0.6125312392 

17.689447 

10.835339 

48.197989 

0.081 

6 

0.1047198 

0.61290441787 

31 .447906 

19.274561 

85.737517 

0.081 

6 

0 .1047198 

0.61338494932 

49.137353 

30.140113 

134.0699 

0.081 

6 

0 .1047198 

0 .61205216674 

0 

0 

0 

0.081 

8 

0. 1396263 

0.78914047544 

1 .9654941 

1 .551051 

6.8994185 

0.081 

8 

0 . 1396263 

0.78935331564 

7.8619765 

6.2058773 

27.6051 17 

0.081 

8 

0.1396263 

0.78970831893 

17.689447 

13.969504 

62.139448 

0.081 

8 

0.1396263 

0.79020589049 

31 .447906 

24.850321 

110.53973 

0.081 

8 

0 . 1396263 

0.79084659909 

49.137353 

38.860109 

172.85838 

0.881 

10 

0.1745329 

0 .9661755943 

1 .9654941 

1 .8990125 

8.4472283 

0.081 

10 

0 .1745329 

0 .96644164455 

7.8619765 

7.5981415 

33.798217 

- 0 .081 

10 

0 .1745329 

0 .96688539866 

17.689447 

17. 103668 

76.080907 

0.081 

10 

0.1745329 

0 .96750736312 

31 .447906 

30.426081 

135.34195 

0.081 

10 

0 .1745329 

0.96830824887 

49.137353 

47.580105 

21 1 .64685 



riiiif ip 


Re wino 


ReUT 


Ue i oh t 
guess (N) 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
55 
- 55 

55 
55 
55 
55 


Ki nemat i c 

v i scos i ty 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0 .0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 

0.0001572 


C-f skin 
•friction 


1 wi no 1 VT 

<_£t) <±t> 


1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 

1 .46 

0.43175 


0 0 
121862.1 36036.961 
243724.21 72073.922 
365586.31 108110.88 
487448.41 144147.84 
609310.52 180184.81 
0 0 

121862.1 36036.961 
243724.21 72073.922 
365586.31 108110.88 
487448.41 144147.84 

609310.52 180184.81 

0 0 

30465.526 9009.2403 

60931 .052 18018.481 
91396.578 27027.721 

121862.1 36036.961 
152327.63 45046.201 
182793.16 54055.442 
213258.68 63064.682 
243724.21 72073.922 
274189.73 81083.163 
304655.26 90092.403 
335120.79 99101 .643 
365586.31 108110.88 
396051 .84 117120.12 
426517.36 126129.36 
456982.89 135138.6 

487448.41 144147.84 
517913.94 153157.08 
548379.47 162166.33 
578844.99 171175.57 

609310.52 180184.81 

0 0 

121862.1 36036.961 
243724.21 72073.922 
365586.31 108110.88 
487448.41 144147.84 
609310.52 180184.81 
0 0 

121862.1 36036.961 
243724.21 72073.922 
365586.31 108110.88 
487448.41 144147.84 
609310.52 180184.81 
121862.1 36036.961 
243724.21 72073.922 
365586.31 108110.88 
487448.41 144147.84 
609310.52 180184.81 


Rg nacel le Cf wing C-f VT C± nacel JeFF -factors 

0 ERR ERR ERR 

121862.1 0.0068483 0.0090923 0.0068483 

243724.21 0.0059038 0.0077087 0.0059038 

365586.31 0.0054333 0.0070316 0.0054333 

487448.41 0.0051303 0.0066 0.0051303 

609310.52 0.0049109 0.0062899 0.0049109 
0 ERR ERR ERR 

121862.1 0.0068483 0.0090923 0.0068483 

243724.21 0.0059038 0.0077087 0.0059038 

365586.31 0.0054333 0.0070316 0.0054333 

487448.41 0.0051303 0.0066 0.0051303 

609310.52 0.0049109 0.0062899 0.0049109 

0 ERR ERR ERR 

30465.526 0.009479 0.0131057 0.009479 

60931.052 0.0080159 0.0108458 0.0080159 
91396.578 0.0073022 0.0097685 0.0073022 

121862.1 0.0068483 0.0090923 0.0068483 
152327.63 0.0065228 0.0086118 0.0065228 
182793.16 0.0062725 0.0082449 0.0062725 
213258.68 0.0060711 0.0079514 0.0060711 

243724.21 0.0059038 0.0077087 0.0059038 
274189.73 0.0057615 0.007503 0.0057615 

304655.26 0.0056382 0.0073254 0.0056382 
335120.79 0.0055298 0.0071698 0.0055298. 

365586.31 0.0054333 0.0070316 0.0054333 
396051.84 0.0053466 0.0069077 0.0053466 
426517.36 0.005268 0.0067957 0.005268, 

456982.89 0.0051962 0.0066936 0.0051962 

487448.41 0.0051303 0.0066 0.0051303 

517913.94 0.0050694 0.0065137 0.0050694 
548379.47 0.0050128 0.0064338 0.0050128 
578844.99 0.0049602 0.0063594 0.0049602 

609310.52 0.0049109 0.0062899 0.0049109 

0 ERR ERR ERR 

121862.1 0.0068483 0.0090923 0.0068483 

243724.21 0.0059038 0.0077087 0.0059038 

365586.31 0.0054333 0.0070316 0.0054333 

487448.41 0.0051303 0.0066 0.0051303 

609310.52 0.0049109 0.0062899 0.0049109 
0 ERR ERR ERR 

121862.1 0.0068483 0.0090923 0.0068483 

243724.21 0.0059038 0.0077087 0.0059038 

365586.31 0.0054333 0.0070316 0.0054333 

487448.41 0.0051303 0.0066 0.0051303 

609310.52 0.0049109 0.0062899 0.0049109 

121862.1 0.0068483 0.0090923 0.0068483 

243724.21 0.0059038 0.0077087 0.0059038 

365586.31 0.0054333 0.0070316 0.0054333 

487448.41 0.0051303 0.0066 0.0051303 

609310.52 0.0049109 0.0062899 0.0049109 


t/c wing 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0 .15 

0.15 

0.15 

0 .15 

0.15 

0.15 

0.15 

0.15 

0 .15 

0 .15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0.15 

0 .15 

0.15 

0.15 

0.15 

0.15 

0.15 

0 .15 

0.15 

0 .15 

0.15 

0.15 

0.15 

0.15 

0.15 

0 .15 


t/c VT 


0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 

0.1011 



d nace 1 1 e 

f nacel 1 e 

FF nacel 1 e 

Q wing 

Q VT 

Q nacel 1 e 

Sw wing 

— 

(ft) 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

<f t*2) 
19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

E. . = 

0.42 

3.4761905 

1 . 1006849 

1 

1 

1 .3 

19.413481 

Iw 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

l 

1 .3 

19.413481 

I 1 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

w 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

- - 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

! i 

mJ 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

: ■ 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

13 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

i...; 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

y 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

- - 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

*=» 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

13 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

» 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

= , 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

— 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

fi 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

w- : l 

m 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

E : i 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

i 1 
r i 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

N 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

Li 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 


0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 

- 

0.42 

3.4761905 

1 .1006849 

1 

1 

1 .3 

19.413481 



Sui VT Sw nacel 1 e Cp wino Cp j/T Cp nacelleCp Total/S Leak & 


<f t‘2) 
1.17 

t‘2) 

0.9632123 

ERR 

ERR 

ERR 

ERR 

Protrub 

ERR 

1.17 

0.9632123 

0 .1021241 

0.0151776 

0.0188774 

0 .0141696 

0.0002834 

1.17 

0.9632123 

0.0997389 

0.0145778 

0.016274 

0.0135881 

0.0002718 

1 .17 

0.9632123 

0.0987397 

0.0143042 

0.014977 

0.0133207 

0 .0002664 

1.17 

0.9632123 

0.0981872 

0.0141398 

0.0141416 

0.0131592 

0.0002632 

1.17 

0.9632123 

0.0978405 

0.0140276 

0.0135369 

0.0130486 

0.000261 

1.17 

0.9632123 

ERR 

ERR 

ERR 

ERR 

ERR 

1 .17 

0.9632123 

0 .1021241 

0.0151 776 

0.0188774 

0 .0141696 

0.0002834 

1.17 

0.9632123 

0.0997389 

0.0145778 

0.016274 

0.0135881 

0.0002718 

1.17 

0.9632123 

0.0987397 

0.0143042 

0.014977 

0.0133207 

0 .0002664 

1 .17 

0.9632123 

0.0981872 

0.0141398 

0.0141416 

0.0131592 

0.0002632 

1.17 

0.9632123 

0.0978405 

0.0140276 

0.0135369 

0.0130486 

0.000261 

1 .17 

0.9632123 

ERR 

ERR 

ERR 

ERR 

ERR 

1 .17 

0.9632123 

0.1101375 

0.0170458 

0.0261288 

0.0159523 

0.000319 

1 .17 

0.9632123 

0.1055143 

0.0159809 

0.0220959 

0.0149409 

0.0002988 

1.17 

0.9632123 

0.1033977 

0.0154834 

0.0201286 

0.0144642 

0.0002893 

1.17 

0.9632123 

0 .1021241 

0 .0151776 

0.0188774 

0.0141696 

0.0002834 

1.17 

0.9632123 

0.1012556 

0.0149646 

0.01798 

0.0139637 

0.0002793 

1.17 

0.9632123 

0 .1006187 

0.014805 

0.0172901 

0.0138091 

0 .0002762 

1 .17 

0.9632123 

0.1001288 

0 .0146796 

0.016735 

0.0136873 

0.0002737 

1 .17 

0.9632123 

0.0997389 

0.0145778 

0.016274 

0.0135881 

0.0002718 

1.17 

0.9632123 

0.0994206 

0.0144929 

0.0158817 

0.0135054 

0.0002701 

1 .17 

0.9632123 

0.0991557 

0.0144208 

0.0155418 

0.0134349 

0.0002687 

1.17 

0.9632123 

0.0989316 

0.0143586 

0.015243 

0.0133741 

0.0002675 

1 .17 

0.9632123 

0.0987397 

0.0143042 

0.014977 

0.0133207 

0 .0002664 

1.17 

0.9632123 

0.0985735 

0.0142561 

0.0147379 

0.0132736 

0.0002655 

1 .17 

0.9632123 

0.0984283 

0.0142132 

0.0145212 

0.0132315 

0 .0002646 

1.17 

0.9632123 

0.0983005 

0.0141746 

0.0143233 

0.0131936 

0.0002639 

1 .17 

0.9632123 

0.0981872 

0.0141398 

0.0141416 

0.0131592 

0.0002632 

1 .17 

0.9632123 

0.0980862 

0.014108 

0.0139737 

0.013128 

0 .0002626 

1 .17 

0.9632123 

0.0979957 

0.0140789 

0.0138179 

0.0130993 

0.000262 

1.17 

0.9632123 

0.0979142 

0.0140522 

0.0136727 

0.0130729 

0.0002615 

1 .17 

0.9632123 

0.0978405 

0.0140276 

0.0135369 

0.0130486 

0.000261 

1 .17 

0.9632123 

ERR 

ERR 

ERR 

ERR 

ERR 

1 .17 

0.9632123 

0 .1021241 

0.0151776 

0.0188774 

0 .0141696 

0.0002834 

1 .17 

0.9632123 

0.0997389 

0.0145778 

0.016274 

0.0135881 

0.0002718 

1 .17 

0.9632123 

0.0987397 

0.0143042 

0.014977 

0.0133207 

0 .0002664 

1 .17 

0.9632123 

0.0981872 

0.0141398 

0.0141416 

0.0131592 

0.0002632 

1.17 

0.9632123 

0.0978405 

0 .0140276 

0.0135369 

0.0130486 

0.000261 

1.17 

0.9632123 

ERR 

ERR 

ERR 

ERR 

ERR 

1 .17 

0.9632123 

0 .1021241 

0.0151776 

0.0188774 

0 .0141696 

0.0002834 

1.17 

0.9632123 

0.0997389 

0.0145778 

0.016274 

0.0135881 

0.0002718 

1.17 

0.9632123 

0.0987397 

0.0143042 

0.014977 

0.0133207 

0 .0002664 

1 .17 

0.9632123 

0.0981872 

0.0141398 

0.0141416 

0.0131592 

0.0002632 

1 .17 

0.9632123 

0.0978405 

0.0140276 

0.0135369 

0.0130486 

0.000261 

1.17 

0.9632123 

0 .1021241 

0.0151776 

0.0188774 

0 .0141696 

0.0002834 

1 .17 

0.9632123 

0.0997389 

0.0145778 

0.016274 

0.0135881 

0.0002718 

1.17 

0.9632123 

0.0987397 

0.0143042 

0.014977 

0.0133207 

0 .0002664 

1 .17 

0.9632123 

0.0981872 

0.0141398 

0.0141416 

0.0131592 

0.0002632 

1.17 

0.9632123 

0.0978405 

0.0140276 

0.0135369 

0.0130486 

0.000261 



CDo 

£ 

K 

CD i 

**CD** 

**** 


- 

cru i se 

cruise 

ERR 

1 

0.0454728 

0.0002983 

ERR 

0.014453 

1 

0.0454728 

0.0002983 

0.0147514 

0.013859? 

1 

0.0454728 

0.0002983 

0.0141583 

0.0135872 

1 

0.0454728 

0 .0002983 

0.0138855 

0.0134224 

1 

0.0454728 

0.0002983 

0.0137208 

0.0133095 

1 

0.0454728 

0.0002983 

0.013607? 

ERR 

1 

0.0454728 

0 .0030273 

ERR 

0.014453 

1 

0.0454728 

0.0030277 

0.0174807 

0.013859? 

1 

0.0454728 

0.003028? 

0.0168888 

0.0135872 

1 

0.0454728 

0.003031 

0 .0166182 

0.0134224 

1 

0.0454728 

0.003033? 

0 .0164564 

0.0133095 

1 

0.0454728 

0.0030377 

0.0163472 

ERR 

1 

0.0454728 

0 .008606 

ERR 

0.0162714 

1 

0.0454728 

0.0086061 

0.0248774 

0.0152397 

1 

0.0454728 

0 .0086063 

0.023846 

0.0147534 

1 

0.0454728 

0 .0086068 

0.0233602 

0.014453 

1 

0.0454728 

0.0086074 

0.0230604 

0.014243 

1 

0.0454728 

0 .0086082 

0.0228512 

0.0140852 

1 

0.0454728 

0.0086091 

0.0226944 

0.013961 

1 

0.0454728 

0.0086103 

0.0225713 

0.013859? 

1 

0.0454728 

0 .00861 16 

0.0224715 

0.0137755 

1 

0.0454728 

0.0086131 

0.0223886 

0.0137036 

1 

0.0454728 

0.0086147 

0.0223184 

0.0136415 

1 

0.0454728 

0 .0086166 

0.0222581 

0.0135872 

1 

0.0454728 

0 .0086186 

0.0222058 

0.01353? 

1 

0.0454728 

0.0086208 

0.022159? 

0.0134961 

1 

0.0454728 

0.0086232 

0.0221193 

0.0134574 

1 

0.0454728 

0.0086257 

0.0220832 

0.0134224 

1 

0.0454728 

0.0086285 

0.022050? 

0.0133905 

1 

0.0454728 

0.0086314 

0.022021? 

0.0133613 

1 

0.0454728 

0.0086345 

0.0219957 

0.0133344 

1 

0.0454728 

0.0086377 

0.0219721 

0.0133095 

1 

0.0454728 

0.0086412 

0.0219507 

ERR 

1 

0.0454728 

0.0170345 

ERR 

0.014453 

1 

0.0454728 

0.0170374 

0.0314905 

0.013859? 

1 

0.0454728 

0.0170463 

0.0309062 

0.0135872 

1 

0.0454728 

0.0170612 

0.0306483 

0.0134224 

1 

0.0454728 

0.017082 

0.0305044 

0.0133095 

1 

0.0454728 

0.0171088 

0.0304183 

ERR 

1 

0.0454728 

0.0170345 

ERR 

0.014453 

1 

0.0454728 

0.028317? 

0.042770? 

0.0138599 

1 

0.0454728 

0.0283332 

0.0421931 

0.0135872 

1 

0.0454728 

0.0283586 

0.0419458 

0.0134224 

1 

0.0454728 

0.0283944 

0.0418168 

0.0133095 

1 

0.0454728 

0.0284405 

0.04175 

0.014453 

1 

0.0454728 

0.0424487 

0.0569017 

0.0138599 

1 

0.0454728 

0.0424721 

0.056332 

0.0135872 

1 

0.0454728 

0.0425111 

0.0560982 

0.0134224 

1 

0.0454728 

0.0425658 

0.0559882 

0.0133095 

1 

0.0454728 

0 .0426363 

0.0559458 



— L/D — 


Drag cru 
(Ibf) 
ERR 

0.0023485 
0.0090163 
0.0198958 
0 .0349506 
0.054161 1 
ERR 

0.0088656 

0.0342689 

0.0758953 

0.1336756 

0.2076117 

ERR 

0.0013295 
0.0050975 
0.0112361 
0.0197196 
0.0305337 
0.0436693 
0.0591203 
0.0768828 
0.0969541 
0.1193328 
0.1440183 
0.1710105 
0.2003104 
0.2319191 
0.2658385 
0.3020707 
0.3406184 
0.3814846 
0.4246727 
0 .4701864 
ERR 

0.0378858 

0.1487707 

0.3320849 

0.5879584 

0.9168104 

ERR 

0.0663398 
0.261845 
0.5859621 
1 .0391612 
1 .622409 
0.108057 
0.4280183 
0.9594856 
1 .7035016 
2.6619076 


Praia 

<N> 

ERR 

0 .0104466 
0.0401061 
0.0885004 
0 .1554674 
0.2409192 
ERR 

0.0394361 

0.1524348 

0.3375977 

0.5946158 

0.9234985 

ERR 

0.0059138 
0 .0226748 
0.0499803 
0.0877167 
0.1358199 
0 .1942497 
0.262979 
0.3419901 
0.4312714 
0.5308163 
0 .640622 
0.7606891 
0.8910208 
1 .0316227 
1 .1825028 
1 .3436709 
1 .5151387 
1 .6969197 
1 .8890289 
2.0914831 
ERR 

0 .1685238 
0 .6617616 
1 .4771803 
2.6153566 
4.0781559 
ERR 

0.2950928 
1 .1647387 
2.6064767 
4.6223969 
7.2167996 
0.4806593 

I . 9039108 
4.2679839 
7.5775158 

I I . 840697 


cruise 

ERR 

5.4910191 

5.7210455 

5.8334202 

5.9034604 

5.9524358 

ERR 

14.761151 

15.281598 

15.53583 

15.696157 

15.810737 

ERR 

17.487202 

18.243889 

18.623756 

18.86656 

19.040184 

19.172799 

19.278614 

19.365721 

19.439131 

19.502145 

19.557039 

19.605447 

19.648576 

19.687337 

19.72244 

19.754441 

19.783787 

19.810839 

19.835897 

19.859206 

ERR 

19.437805 

19.810403 

19.985801 

20.092345 

20.165012 

ERR 

18.45041 

18.708131 

18.826871 

18.896846 

18.94244 

16.979733 

17.156185 

17.235576 

17.280557 

17.307967 



APPENDIX E 


MATHCAD RESULTS FOR 
STRUCTURAL ANAYLSIS 


VI 


JiliiJii.Lju 


Read Pints from data file (hp33.pm-section points, hp33int.pm-lift, drag, moment) 


i:=0..45 j :=0.. 12 

point : = READPRN( points ) LDM = READPRN( 1dm) 

x u. := P° int i- h 45 ,o *□. : -point i+45 I :=point 45 . x 

*I := *u 


44,0 


:= 0.999 


Lift:=65.769 Span: s 1,25 



zl : - 0 , 0.25 .. 1 .25 N = newton 
Equation for the elliptical distribution of lift across the span of the wing 



Polynomial fit of: Force. Shear and Moment curves 


Rows -9 k-0..Rows Zf :-k-0.125 yj- := y^ZfJ degree :=5 i :=0.. degree point :=0.. Rows 


' matrix 


(point, i) 


: f. 


_ ( t y 1 

\ z matrix ' z matrix/ 1 


(point, 0)J y fjcoeff \ z matrix z matrix/ ' z matrix 


P 


y f_eq( z ) ’ I y f_coeff. z 


A 


N 


Force 



z fk- zl 






Vf :=v^ZfJ degree :=5 i: = 0..degree point =0.. Rows 
matr ’ X (poi“*.') . ^(point.o) v f_coeff " ( z matrix z matrixj ' z matrix -v f 

v f_eq^ 2 ^ ~ ! v f.coeff.- z |- N 


T - v 1 .- T 

rix/ 


Shear 



v f coeff - 


-0.003 

5.102 

42.577 

-44.308 

28.794 

-8.037 


z f. ,zl 


M 


v =M M 


degree : = 5 i : = 0 .. degree point : = 0 .. Rows 


z matrix 


(point 


:=fz i 

int.i) 


( point, 0) 

^ f eq^ z ^ f coeff . ' z \ 


M 


._( t y 1 t M 

f_coeff \ z matrix z matrix/ * z matrix f 


/ i 


m 


Moment 



0 


1 




0 


0 

0.00182 


0.01182 

0.00455 


0.0185 

0.01136 


0.0279 

0.02273 


0.03786 

0.03409 

II 

0.04522 

0.04545 


0.05091 

0.06818 


0.05965 

0.09091 


0.06645 

0.11364 


0.0719 

0.13636 


0.07648 

0.15909 


0.08044 

0.18182 


0.08372 

0.22727 


0.08889 

0.27273 


0.09245 

0.31818 


0.09463 

0.36364 


0.09545 

0.42209 


0.09505 

0.46755 


0.09335 

0.513 


0.09015 

0.53573 


0.0877 

0.55845 


0.08522 

0.58118 


0.08194 

0.60391 


0.07822 

0.62664 


0.07404 

0.64936 


0.06945 

0.67209 


0.06451 

0.69482 


0.05925 

0.71755 


0.05375 

0.74027 


0.04745 

0.763 


0.04222 

0.78573 


0.03621 

0.80845 


0.03012 

0.83118 


0.02399 

0.85391 


0.01875 

0.87664 


0.01363 

0.89936 


0.00951 

0.91559 


0.007 

0.93182 


0.005 

0.94319 


0.00425 

0.95455 


0.00375 

0.96591 


0.0033 

0.97727 


0.00415 

0.98864 


0.0055 

r\ r\r\n> 


n r\rvn 


' 0 


0 

- 0.0086 


0.00182 

- 0.01275 


0.00455 

- 0.01812 


0.01136 

- 0.02377 


0.02273 

- 0.02818 

x ,= 

0.03409 

- 0.03152 


0.04545 

- 0.03717 


0.06818 

- 0.04151 


0.09091 

- 0.04493 


0.11364 

- 0.04792 


0.13636 

- 0.05034 


0.15909 

- 0.05134 


0.18182 

- 0.05234 


0.22727 

- 0.05234 


0.27273 

- 0.05234 


0.31818 

- 0.05234 


0.36364 

- 0.05234 


0.42209 

- 0.05234 


0.46755 

- 0.05234 


0.513 

- 0.05134 


0.53573 

- 0.05019 


0.55845 

- 0.04766 


0.58118 

- 0.0447 


0.60391 

- 0.04132 


0.62664 

- 0.03767 


0.64936 

- 0.03383 


0.67209 

- 0.02988 


0.69482 

- 0.02588 


0.71755 

- 0.02174 


0.74027 

- 0.01793 


0.763 

- 0.01402 


0.78573 

- 0.01052 


0.80845 

- 0.00752 


0.83118 

- 0.0055 


0.85391 

- 0.00466 


0.87664 

- 0.00466 


0.89936 

- 0.00475 


0.91559 

- 0.0042 


0.93182 

- 0.004 


0.94319 

- 0.0038 


0.95455 

- 0.003 


0.96591 

- 0.001 


0.97727 

0.001 


0.98864 

/> rv/\f> 


n onn 


Airfoil Coordinates 
Of the HPRS-3A 



\J.7 77 


t 




\J.7 7 7 

1 


0.005 


0.005 


1 


Geometric Parameters 

i:=0..45 j : = 0..44 KPa=1000-Pa 


Span 

s := 1.25-m 

Chord Lenght 

cord :=0.455-m 

Engine Weight 

w e :=5.719-N 

Cross Sectional 

A e : = (0.0455-m) : 

Area of Engine 

A e = 0.007* m 2 

Distance of 


Engine from 


the root chord 

d :=57-cm 

Center of Gravity 

X ecg :=0 - 017 - m 

of the Engine 

Y „„ :=0.006-m 


ecg 

Thrust 

r :=— -N 


30 


Wing Sweep 

P :=0-deg 

Drag of Wing 

D : = 0.109-N 

Moment 
Due to Drag 

M d :=j*(«) M d 

Center of 
Gravity 

X cg :=0.113-m 
Y C g :=0.006-m 

Airfoil 

Coordinates 

Upper 

x a :=x u cord 

i i 

Y u. := yu. cord 

Lower 

Y | :=y j -cord 

Wing Center 
of Gravity 

X weg := 0-209-m 
Y wcg := 000561 m 


Calculated Cross 
Sectional Area 
at Root 

Moment of inertia 
about the x-axis 
upper 

Moment of inetrtia 
about the x-axis 
lower 

Moment of inertia 
about the y-axis 
upper 

Moment of inertia 
about the y-axis 
lower 

Total moment of 
inertia of the wing 
about x-axis 




Total moment of 
inertia of the wing 
about y-axis 


*wy *wyl + *wy2 ^^ws'(®* m ^ wcg) 


Total moment of 
inertia about x-axis 

Total moment of 
inertia about x-axis 


* X ■ ~ * WX ^ 

^ y ~ ^ wy ^ 


ws'(^ 

ws'(^ 



+- A 


+■ A 


e ( Y eg Y ecg) I 

e‘( X cg -X ecg) 1 


•N 


A ws =0.02*m 2 

wxl = 1-772*10 5 *m 4 

wx2 = -3.314* 10~ 6 *m 4 
wyl - 5.537* 10~ 4 *m 4 

wy2 =-2.836*10“ 4 *m 4 

wx =2.168* 10~ 5 *m 4 
wy = 0.002*m 4 

x = 2.168* 10 -5 *m 4 
= 0.002* m 4 


ORIGINAL PAGE IS 
OF POOR QUALITY 


(T-sin(p)) 










o 


KPa 


'3.892 ' 


3.892 

6.397 


0.449 

7.605 


-1.498 

9.168 


-4.237 

10.812 


-7.141 

12.094 

•KPa a 2 = 

-9.287 * 

13.064 


-10.947 

14.705 


-13.497 

15.965 


-15.483 

16.957 


-17.075 

17.823 


-18.414 

18.523 


-19.572 

18.81 


-20.532 

19.092 


-22.047 

19.082 


-23.094 

19.073 


-23.738 

19.064 


-23.986 

19.052 


-23.882 

19.042 


-23.396 

19.033 


-22.473 

18.737 


-21.764 

18.397 


-21.046 

17.656 


-20.096 

16.789 


-19.017 

15.799 


-17.804 

14.732 


-16.471 

13.608 


-15.037 

12.453 


-13.51 

11.283 


-11.912 

10.073 


-10.082 

8.958 


-8.563 

7.815 


-6.817 

6.79 


-5.048 

5.912 


-3.267 

5.319 


-1.745 

5.069 


-0.258 

5.065 


0.937 

5.088 


1.665 

4.924 


2.244 

4.863 


2.46 

4.803 


2.604 

4.567 


2.732 

3.983 


2.482 

3.398 


2.087 

2.813 


1.648 

2.23 


2.23 


Deflection Angle 


a0 ■ - M fjcoeff^'N-m 

a3:=M f_coeff 3 ' N - m 

M we :rw e d 

al - M fjcoeff^N^m 

a4 -Mf coeff 

“ 4 


a2 -M ^coeff^N-m 

a5: = M f_coeff s - N - m 



5 4 3 2 

Moment Equation: Mo(x):-a5*x +a4*x +■ a3*x -{-a2-x +- al-x + aO -h M we 

01(x)ss-!-*al-x 2 -{- aOx-K -*a3*x 4 -t- -*a2-x 3 +- i-a5-x 6 + -*a4*x 5 - M wp *x+- Cl 
2 4 3 6 5 e 

01( x) *al*x 2 -h aOx -h -a3x 4 +- -*a2*x 3 + i*a5-x 6 +- — -a4-x 5 - M wp -x 
2 4 3 6 5 we 

y(x)*-!— a5*x 7 + — -a2-x 4 +- — -al-x 3 -h — *a4*x 6 + — *a3-x 5 -h— -aO-x 2 - wp *x 2 + Cl-xi- C2 
42 12 6 30 20 2 2 we 

Cl :=-0l|— j Cl =-2.231 *m*N 

yl(x) :=— -a5-x 7 -i- — -a2-x 4 +- i-al-x 3 +■ — -a4-x 6 -i- — -a3-x 5 + — aO-x 2 -t- Cl-x- — M ,,x 2 
42 12 6 30 20 2 2 we 

C2 :=- yl j— j C2 =3.551*m *N 


y(x) 


:=— -a5-x 7 -t- — -a2-x 4 +- — -al-x 3 +■ — -a4-x 6 -i- — -a3-x 5 -i- — -aO-x 2 - — M ^..-x 2 + Cl x+ C2 

42 12 6 30 20-2 2 we 


X 


- 0 . 01..0 

m m 


Dihedral* 



♦Mesurment of Dihedral is not devided by E, 
the modulus of elasticity of the carbon fiber 
and I, the moment of inertia of the wing 


- 1 . 11 * 100.25 


0.5 


0.75 


1 


1.25 



Shear Flow 


Shear flow analysis modeled as a closed loop section. 

Assuming counterclockwise direction is positive 

carbon = 0.05- 10* 2 *m w e =5.719 *N x C p'=0.111-m Lift = Lift-N 

Y 

w e' A ecg 9 ' x cp , - 

x := x =-1.837* 10 s *m z *N 

2.A ws -carbon ±ick 



APPENDIX F 


INVISID VORTEX PANEL CODE 


VII 


Angle of Attack : 


4 


Flap Angle: -20 


COEFF : Ith panel no. 90 



alpha=4 cl=2. 16848 cd=0 . 00768639 cnle— 9.2478 cnc4=-8 .70568 

I 

Press any key to exit 


This program uses an inviscid method to solve for the 
velocity and pressure coefficient distribution over an airfoil 
shape. The program requires a file which contains the airfoil 
coordinates starting at the trailing edge on the lower side and 
then proceeding around the airfoil all the way back to the to the 
starting point. The second file is angle.dat which contains the 
angles of attack for the airfoil. If flap deflection is required 
a file with the flap deflection angles is needed. The program 
produces two files cp.out, and vr.int. The first contains the 
coefficient of pressure distribution over the airfoil and the 
second contains the lift drag and moment coefficient. 





/*****«*•*********************************************** 
♦••••••♦•a********* ************************* a*********** 

mmmmmum Program Vortex.c ********************** 
************** December 1992 ********************** 
************** Yuri Zhovnirovsky, Shilpa Shroff ******** 


^include <stdio.h> 
jfinclude <malloc.h> 

^include < stdlib .h > 
include <math.h> 

^include <conio.h> 

#include < graphics. h> 

include "c:\mqp\mv.c" /* Numerical Recipes Dynamic Memory Allocation */ 
^include "c:\mqp\inverse.c" /* Calculates an inverse of a matrix */ 

#define FALSE 0 
^define TRUE 1 

float PI; 

int GraphDHver; /* graphics variables for the initialization of */ 

int GraphMode; /* graphics screen */ 

/* This procedure solves a system of linear equations 
c * Contains the variable coefficients (left hand side) 
a - Right hand side 
x - Solution vector 
n - Number of equations 

*/ 

void solve(c, a, x, n) 
float **c, *a, *x; 
int n; 

{ 

float **cc; 

register int i,j,k; 

cc = matrix(l ,n,l ,n); 
inverse(c,cc,n,n); 
mproduct(cc,a,x,n,n,n, I); 
free_matrix(cc, l ,n,l ,n); 

} 

/* This procedure integrated the pressure coefficient in order 
to determine the Lift, Drag and Moment about the quater chord 
point. 

fp - file pointer for output. 

x,y - panel mid point coordinates 

xb, yb - panel end point coordinates 

s * panel lenght 

cp - pressure coefficient 

th - angle of panel 

M - number of panels 

mid_point - Separation point between upper and lower surface 
alpha - angle of attack 

*/ 

void integ(fp,x,y,xb,yb,s,cp,th,M,mid_poini, alpha) 

FILE *fp; 

float *x, *y, *xb, *yb, *s, *cp, *th; 
int M, mid_point; 

float alpha; 

{ 

register int j,k; 

int M2; 

float cn, ca, cmle, cmc4, cnl, cnu, cal, cau, cl, cd; 

float q_inf= 137.99, p_inf= 101325.0; 



cau=0; cal = 0; 

cnu = 0; cnl = 0; 

cmle = 0; 

M2 = mid_point; 

for (j = 1 ; j < M2; j + + ) { 

cnI = cnl + cp[jj*s[jj*cos(th[jl); 
cal = cal + cp|j]*s[]l*stn(th[j|); 
cmle = cmle-(p_inf+cp[j]*q_inf) * s[j] / q_inf * 

(x[jl*cos(PI-th(j]) + fabs(y[j]) * sin{Pl-th[j])); 

} 


for (j = M2; j < =M; j+ +)( 

cnu = cnu + cp{j]*s(j]*cos(th|jl); 
cau = cau + cp[j]*s[j]*sin(th[j]); 
cmle = cmle + (p_inf+cp{j]*q_inf) * sjj ] / q_inf * 

(x{j]*cos(th[j|) + fabs(y[jl) * sin(th[j])); 

} 

printf("cnu = %g cnl= %g cau = %g cal = %g \n",cnu,cni,cau,cal); 

cn=-cnu-cnl; 

ca = cau + cal; 

cl = cn*cos(alpha)-ca*sin(alpha); 
cd — cn*sin(alpha) + ca*cos(alpha); 
cmc4 = cmle + cl/4 .0; 

printf("alpha = %g cl = %g cd = %g cmle = %g cmc4=%g Nn* 1 , 

alpha* 1 80/PI, cl, cd, cmle, cmc4); 

fprintf(fp,"%g, %g, %g, %g, %g\n",alpha*l 80/PI, cl, cd, cmle, cmc4); 

} 

I* Scans the input file to determine the number of panels. 

Also determines the position of the inflection point between 

the upper and lower surface 

fp - File pointer 

mid * mid point (return vector) 

Procedure returns the number of panels. 

*f 

int scan_file(fp, mid) 

FILE * *fp; 
int *mid; 

{ 

int i; 

float xa, xa_old, ya; 

i=l; 

*mid = -l; 

while((fscanf(fp,"%f %f ,&xa, &ya))! = EOF){ 
if (i-=»l) 

xa_old = xa; 

if (((xa_old-xa)<0.0) Sl& (*mid<0)) 

*mid = 
xa_old = xa; 
i+ + ; 

} 

fclose(fp); 

retum(i-l); 

} 

/* Read in the panels from the file 
fp - File pointer 
x_p - x coordinate 
y_p - y coordinate 

*/ 

void getpanels(fp,x_p,y_p) 

FILE ♦Ip; 

float *x_p, *yj>; 



{ 

int i; 

float xa, ya; 


i*= 1 ; 

while(((fscanf(fp," %f %r,&xa,<&ya))!= : EOF)){ 
x_p[i| = xa; 

y_p[i]=ya; 
i + + ; 

} 

fclose(fp); 

} 

/* Rotate preforms a 2-D rotation of the airfoil coordinates 
about the a certain percent of the chord. It only rotates the 
points which are after this percentage. This is done to allow 
for flap deflection analysts, 
r - Degrees of rotation 

x_in, y in - input x, y coordinate vector 
x_out t y_out - rotated x, y coordinate output vector 
p_cbord - percent of chord to rotate about 
Mp - Number of panels 

*/ 

void rotate(r, xjn, yjn, x_out, y_out, p_chord, Mp) 
float r, *xjn, *yjn, *x_out, *y_out, p_chord; 

int Mp; 

{ 

float z_matrix[3](3]; 

float xu_o, yu_o, xl_o, yl_o; 

int t, first_pass; 

/* Z matrix 

•* 

** j cos A -sin A 0 | 

** | sin A cos A 0 | 

** | 0 0 l | 

** ___ 

*/ 


z_mat ri x (0 ] [ 0 1 = cos ( r) ; 
z_matrix[0J[ l J = -sin(r); 

z_matrix[ 1 j[0| = sin(r); 
z_matrix[ 1111] = cos(r); 

i= l; 

xu_o = -10; yu_o = -10; 

xl_o = -10; yl_o = -10; 

p_chord = l-p_chord; 

while ((i < Mp+ I) && (yu_o= =-10)){ 

if ((xjn[i]< =p_chord) &&. (xl_o* =-10)){ 
xl_o = x_infi-l]; 
yt o = y in[i-l]; 

} 

if ((xjnfi] > =p_chord) && (xJ_o! = -lO)){ 
xu_o = x_in[il; 
yu_o = y_m[i); 


i + +; 

} 

for (i = 1; i< = Mp + 1; i+ 4-){ 
if (xjnfi] > =p_chord){ 
’if (i < Mp/2){ 



x_out[i] = xl_o + (x_in[i|-xl_o) * z_matrix[0|[0] + 
(y_in(i]> * z_matrixIO][l ]; 
y_outfil = (x_infi]-x!_o) * z_matrix( l jfO| + 

(y in|i]) * z_m*trixfl](l ]; 

} 

else{ 

x_out(i] = xu_o + (x_in{i|-xu_o) * z_matrixj0)(0] 4- 
(y_infil) * z_matrixf01[ 1]; 
y_out[il = (xjn[i|-xu_o) * z_matrix[ll[0] + 
(y_infi|) ♦ z_matrix[ 1 1[ I]; 

} 

} 

else{ 

x_out(i] = x_in(i]; 

y_outfi] =yjn[ij; 


5 

} 

/* Create Graphical output */ 

void do_graphics(M,xb,yb,x,y) 
int M; 

float *xb, *yb, *x, *y; 

{ 

int i; 

float y_f, x_f, x_o, y_o; 

y_f=(float)getmaxyO/l .5; 

x_f=(float)getmaxxO/l *5; 

x_o = x_f/3.5; 

y_° = y_f/ 1 .4; 

setcolor(RED); 

line(0, y_o, getmaxxO. y_o); 

Une(getmaxx0/4, 0, getmaxxO^*, getmaxyO); 
setcolor(GREEN); 
for (i = 2; i<=M+l; i++){ 
setcolor(BLUE); 

setflllstyIe(SOLrD_FILL, YELLOW); 
if (i < =M) 

filieiIipse(x[i]*x_f+x_o, -y[il*y_f+y_°> 2, 2); 
setcolor(GREEN); 

line(xbfi-l]*x_f+x_o,-yb[i-l|*y f+y o,xb[i]*x f+x_o,-yb[i]*y_f+y_o); 

} 

line(xb[l|*x f+x_o»-yb[l]*y_f+y o,xb[M+II*x f+x o,-yb[M+ 1 ]*y_f+y_o); 

} 


void main0{ 


register int 


int 

M, mpl, ipl, mid_point; 

short int 

do_flaps; 

float 

angle, flap_angle, p_ch; 

FILE 

*fp, fp_out, *fp_int, *fp_angle; 

char 

filename_air[I5J, fllenamef 15}, ch; 

float 

A, B, C, D, E, F, G, P, Q, alpha; 

float 

*xb, *yb, *xbl, *ybl, *x, *y, *$, 


♦theta, *v, *cp, ♦rhs, *garru 
♦♦cn2, ♦*ctl, **ct2, ♦♦an, ' 


'sine, ^cosine, 


clrscrO; 

do_flaps = FALSE; 


printf("Enter input file name: "); 
scanf( " % s " , filena mej ir); 



fflush(stdin); 

fp = fopen(fUename_air, "r"); 
if (fp= = NULL)} 

printf("FILE NOT FOUNDSn”); 
exit(-l); 

j 

printf(*Do you wish to deflect the flap: [y/n]"); 
ch=getch(); 

if((ch*»’y’) || (ch= = , Y’)){ 

printf("\n Enter the file name with the flap deflection: *); 

scanf(" %s", filename); 

fflush(stdin); 

fp_angle = fopen(filename, *r*); 
if (fp_angle= = NULL){ 

printf("ERROR: Flap file not found\n"); 
exit(-l); 



} 

printf("\n Enter the angle of attack: "); 
scanf(" % r ,<fcf1ap jingle); 
fflush(stdin); 

printf("\n Enter the size of the flap in percent chord: "); 

scanf(’ , %r,&p_ch); 

fflush(stdin); 

p_ch = p_ch/ 100.0; 

do flaps = TRUE; 

} 

else 

printfC\n"); 

M — scan_file(fp, &mid_point); 
fp = fopen(filename jiir, V*); 

I* Allocate memory for coordinates */ 

xbl =vector(l ,M); 
ybl =vector(l,M); 
xb = vector(l,M); 
yb=vector(l f M); 

M-; 

printf("Using %d paneIs\n w ,M); 

getpanels(fp,xbl ,ybl); /* Read in the panels */ 

/• Allocate memory */ 

x=vector(l ,M); 
y=vector(l ,M); 
s=vector(I ,M); 
sine = vector(I ,M); 
cosine =vector( 1 ,M); 
theu = vector( 1 ,M); 
rhs=vector(l ,mpl); 
pi = 4.0*atan(1.0); 
mp l = M + l ; 

sprintf(filename, "vr.int"); 

fp_int = fopen(filename,"w"); /* open output file vr.int 

This file will contain the 
integrated cp values */ 

fp = fopen("cp.out Vw"); /* Open output file cp.out This will 

contain the raw cp values */ 

if (!do_flaps){ 

fp jingle = fopen( "angle.dat " f "r">; 
if (fp_angle= — NULL){ 

printf(" ERROR: The file (angle.dat) which contains the angles of attack is not found\n"); 
exit(0); 

} 



} 

angle = 0.0; 

/* rotate if needed */ 

route(angle, xbl, ybl, xb, yb, p_ch, M); 

/* Start Graphics */ 

GraphDriver = DETECT; 

initgraph(&GraphDriver, &GraphMode, "c:\\borlandc\\bgi"); 

/* Read in flap defection if needed */ 
whiie(fscanf(fp_angle,"%P\&angle)! = EOF){ 
alpha = angle* PI/ 1 80 .0; 
if (do_flaps){ 

gotoxy(30,4); printf("F1ap Angle: %g ".angle); 
rotate(alpha, xbl, ybl, xb, yb, p_ch, M); 
alpha = flap_angle*PI/l 80.0; 
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gotoxy(5,4); printfC Angle of Attack: %g ".alpha* 180.0/PI); 

/* Calculate center points x,y length of panel s 
Angle of panel theta, sine and cosine of this angle, 
and the right hand side matrix. 

*/ 

for (i = l ; i < = M ; i + + ) { 
ip 1 = i + 1 ; 

xfi] =0.5*(xbfij + xbfipl )); 
y[il = 0.5*(yb[ij+ybfipl j); 

sfi] = sqrt(pow(xb[ipl ]-xb[i}, 2) + pow(ybfip 1 1-ybfi], 2)); 

theta(i] = atan2(yb[ipl |-ybfi], xbfip 1 ]-xb(i]); 

sine[i] = sin(theu[ij); 

cosinefi] = cos(theta(i|); 

rhsfij = sin(thetafi)-alpha); 

} 

do _graphics(M,xb,yb,x,y); /* Display graphics *1 

/* Allocate more memory */ 

cnl =matrix(! ,M,1,M); 

cn2 = matrix(l ,M,1 ,M); 

ctl = matrix(l,M, 1 ,M); 

ct2 = matrix(l,M, 1 ,M); 


gotoxy(15,8); printf("COEFF: Ith panel no. %d ",i); 

for(j=l;j<=M;j++){ 

if (i==j){ 

cnl ft] [j ] *-1.0; 
cn2[i)(j]= 1.0; 
ctl[ij(ji = 0.5*PI; 
ct2fi](j|=0.5*PI; 

} 

else{ 

/* Calculate coefficient for the integral *1 


A = -(x[i)-xb[j])*cosine(j]-(yfil-yb{j])*sine(j J; 

B = pow(x[i}-xb(j], 2) + pow(yfi]-yb(j j, 2); 

C = sin(thetafi]-theta(j]); 

D = cos(theu[ij-theta[j]); 

E = (x[i]-xb[ji)*sine[jl-(y[i]-yb[jI)*cosine[j]; 

F = log(l+s[j]*(s[j| + 2.0*A)/B); 

G = atan2(E*s(j], B + A*s(j]); 

P * (xfi]-xb[jl)*sin(theta[i]-2.0*theuij]) + 
(y[ij-yb[j])*cos(theta[il-2.0*theta[jj); 

Q = (x[i]-xb(j])*cos(thetafi)-2.0*theta[j])- 
(y[i]-yb[j])*sin(theufi]-2.0*theta[j]); 
cn2fi][j] = D+0.5*Q*F/s[j]-(A*C + D*E)*G/s[j]; 
cniniU! = 0.5*D*F+C*G-cn2ni(j]; 
ct2fi][j] = C4-0.5*P*F/sljl -f (A*D-C*E)*G/Rfjl; 
ctl fillil =0.5*C*F-D*G-ct2[i)[j]; 



1 

} 

J 

f* Allocate more more more and some more memory *1 
an = matrix(l ,mpl , I ,mp I); 
at = matrix(l,mpl ,l,mpl); 

/* Fill Ld left hand side matrix •/ 

for (i=l; i < = M; i+ +){ 
anfiim-caimm; 
anfijfmpl J = cn2[il[M]; 

®tfij[ 1 1 — ct 1 fil [ 1 1; 
at[i][mpl ]=ct2[iirMl; 
for (j=2; j < = M; j 4- +){ 

an[illj]=cnlfil(j] + cn2fi]0-l]; 

•tmDl«ctl[ilDl+ct2[i]D-ll; 

} 

} 

anfmpl |[1 j = 1 .0; 
anfmpl ][mp l] = 1 .0; 
for (j = 2; j < = M; j + +) 
anfmpl ][j| =0.0; 
rhsfmpl]=0.0; 

/* Free (quickly) some of that allocated memory */ 

free_matrix(cn 1 , 1 ,M, 1 ,M); 

free_matrix(cn2,l ,M, 1 ,M); 

firee_matri x(ct 1 , 1 1 M , 1 , M) ; 

free_matrix(ct2, 1 ,M,1 ,M); 

/* allocate memory for solution vector */ 
gama = vector( 1 , mp l ) ; 

/• SOLVE THE EQUATIONS */ 
solve(an, rhs, gama, mpl); 

free_matrix(an, 1 ,mp 1 , 1 ,mp 1 ); 
v = vector(l ,M); 
cp=vector(l ,M); 

/* Write cp, and velocity to the cp.out file */ 
fprintf( fp, "alpha = %g\n", angle); 
for (i= 1; i < =M; i + +){ 

v[i] = cos(theta[i]-alpha); 
for (j = l ; j < = mp l ; j + + ){ 

v fi] — v[i] 4-at[i][j]*gama(j]; 
cp[i} = 1 ,0-pow(vfi|,2); 

} 

fprintf(fp,"%g, %g, %g, %g, %g, %g\n\ xbfi], ybfil, xfi], yfi], 

thetafi], cpfif); 

} 

fprintf(fp,"\n\n"); 

/* Integrate the Cp coefficients to determine the lift and drag */ 
integ(fp_int, x,y,xb,yb,s,cp, theta, M, mid jjoint, alpha); 

/* FREE FREE FREE the memory */ 

free_matrix(at, 1 ,mpl , I ,mpl); 
free_vector(gama,I ,mpl); 
free_vector(cp, 1 ,M); 
free_vector(v, I ,M); 

} 

/• Close files */ 
fclose(fp); 
fclose(fp_int); 
fclose(fp_angIe); 


f* FREE ALL OF THE TAKEN MEMORY */ 

free_vector(x, 1 ,M“) ; 
free_vector(y, I ,M); 
free_vector( theta, I ,M); 
free_vector(cosine, I ,M); 
free_vector(sine, 1 ,M); 
firee_vector(rhs, I ,mp 1); 
free_vector(s, 1 ,M); 
firee_vector(xb, 1 ,mpl); 
free_vector(yb , 1 , mp 1) ; 
free_vector(xb l , 1 ,mp l); 
free_vector(yb [ , ] ,mp I); 
printf("Press any key to exit\n"); 
getchO; 

cIosegraphQ; /* close graphics */ 


/a***************************** INVERSE .C ************************♦**********/ 

float inverse(in_matrix,out_matrix, rows, columns) 
float **tn_matrix, **out_ matrix; 
int rows, columns; 

{ 

int i,c,r,row, column; 

float temp, factor; 

if (in_matrix[ l]{ 1 1 = = 0){/* check if the first entry is zero */ 

row = 2; 

while (in_matrix|row][l] = = 0){/* find the next non zero entry *f 
row + + ; 

} 

for (cotumn=l; column < = columns; column + + ){ /* switch lines around */ 
temp = in_matrix[row Hcolumn) ; 
in_matrix[rowl [column] — in_matrix[row][column|; 
in_matrix [row] [column] = temp; 



/* create an identty matrix *f 

for (row = I ; row < = rows; row -l- +) 

for (column=i; column < ^columns; column ++) 
if (column = = row) 

out_matrix(row][column] = 1 .0; 

else 

oui_matrixlrow]lcolumn] = 0.0; 

I* since the inverse matrix must be sqare only one counter is used */ 

for (1= I; I < = rows; i + + ){ 

gotoxy(l ,19); printf(" Inverse up: %d \n",i); 

/* introduce leading l’s */ 
if (in_matrix(i][il!- 1){ 

factor= 1.0 / in_matrix[i|[i]; 

for (c= 1 ; c < =columns; c + +){ 

in_matrix[i][c] = in_matrix[i][c]*factor; 
out_matrix(i][cl = out_matrix[il(cl*factor; 

} 

} 

/* intrudce 0 bellow the leading l’s */ 
for (r=i+ 1 ; r < = rows; r+ +){ 
if (in_matrix(r][il! = 0){ 
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factor = in_matrix(rlfil / in_matrix[i|(i|; 

/* subtract the factor */ 

for (c= I ; c< ^columns; c+ + ){ 

in_matrix[r|[cl = in_matrix(rl(c]-factof*in_matrix[i|(c); 
out matrix(r][c|==out_matrix(r][c|-factor*oui matrix[il(cj; 

) 

} 

} 

} 


for (i = rows; i > = l ; i— ) { 

gotoxy(30, 19); print f(" Inverse down: %d \n",i); 

/•introduce 0 above the leading l’s */ 
for (r=i-l; r> = 1; r-){ 

if (in_matrix(r|fil! =0){ 

factor = in_matrix[r][i] / in_matrix(i]fil; 

/• subtract the factor */ 

for (c = columns; c > = 1 ; c— ) { 

in_matrix( r| [c ) = in_matrix[ rj [c ]-factor*in_matrix [i | (c ) ; 
out_matrix(r][c| = out_matrix[rj[c]-factor*out matrix[i|[cj; 

} 

) 

} 


retum(O); 

} 


float mproduct(in_matrixl , in_matrix2 > out__matrix, rows 1 , columns I , rows2, columns2) 
float **in_matrixl, *in_matrix2, *out_matrix; 
int rowsl, columns 1, rows2, columns2; 

{ 

int row l, row2; 

for (rowl = 1; row I < = rowsl; row! + + ){ 
out_matri x [ row 1 ] = 0 ; 

for (row 2 = 1 ; row2 < * rows2; row2 + + ) { 
out_matrix[rowl ] = out_matrixfrowl ] ■+* 
in matrix 1 [row l ][row2] * in matrix2frow2); 

} 

} 

return (0); 


/**************************** fvfV.C *+*++*+*+*+++*********++***/ 

void nrerror(error_text) 
char error text(l; 

{ 

void exitO; 

fprintf(stderr, "Numerical Recipes run-time error. .An"); 
fprintf(stderr, " %s\n" ,error_text); 
fprintf(stderr,"...now exiting to system. .An"); 

*exit(l); 

} 


float *vector(nl,nh) 
int nl.nh; 

{ 



float *v; 


v = (float *)ma Hoc ((unsigned) (nh-nl+ l)*sizeof(float»; 
if (!v) nrerror( "allocation failure in vectorO"); 
return v-nl; 

} 

float **matrix(nrl,nrh,ncl,nch) 
int nrl»nrh.ncl,nch; 

{ 

int i; 

float **m; 

m= (float **) mailoc((unsigned) (nrh-nrl+ l)*sizeof(float*)); 
if (!m) nrerror( "allocation failure I in matrixO"); 
m -= nrl; 


for(i = nrl;i < -nrh;i+ 4-) { 

mfij = (float *) malloc((unsigned) (nch-ncl+ I)*sizeof(float)); 
if (!m[i]) nrerror( "allocation failure 2 in matrixO"); 
m[i] -= ncl; 

} 

return m; 

} 

void free_vector(v,nl,nh) 
float *v; 
int nl,nh; 

{ 

free((ch«r*) (v + nl)); 

} 

void free_matrix(m,nrKnrh,ncl t nch) 

float **m; 

int nrl,nrh,ncl,nch; 

{ 

int i; 

for(i = nrh;i> = nri;i — ) firee((char*) (m[il+ncl)); 
free((char*) (m + nrl)); 

} 
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